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in  the  fieestream  to  fall  into  a  flat  plate  laminar  boundary  layer  and  impact  the  plate  are  alao  discussed.  A  combination  of  dye  flow  visualization  and  a  video 
baaed  particle  tracking  technique  was  used  to  study  the  convecting  particle  problem.  Although  the  Reynolds  number  at  the  critical  conditioo  for  turbulent 
wedge  generation  by  fixed  particles  and  turbulent  spot  gerKration  by  convecting  particles  are  similar,  transition  in  these  two  situationa  appears  to  be 
fundamentally  different.  The  development  of  a  turbulent  wedge  near  the  critical  condition  is  a  relatively  gradual  process.  In  contrast,  turbulent  spou  form 
relatively  quickly  after  the  convecting  particles  enter  the  bountlary  layer  and  impact  the  plate.  Turbulent  wedge  formation  downstream  of  a  fixed  particle 
reaulu  from  the  deaubilization  of  the  near  wall  flow  by  the  vortical  structures  shed  into  particle  wake.  This  shedding  process  is  dominated  by  periodically 
shed  loop  shaped  hairpin  vortices.  Observation  of  subharmonic  oacillatiom  at  Vt  and  U  of  this  shedding  frequency  suggest  that  a  chaotic  route  to  turbulence 
by  a  aeries  of  period  doublii^  bifurcations  is  possible.  Convecting  particle  wakes  have  little,  if  any,  affect  on  turbulent  spot  gerKration.  The  diaturbatKes 
to  the  flow  tKar  the  wall  due  to  the  particle  impact  rapidly  develop  and  spots  begin  to  etiKrge  at  the  particle  in  all  cases. 

The  analysis  first  predicts  the  trajectories  of  the  convecting  particles.  The  development  of  wave  packets  from  disturbances  at  different  locations  is  considered . 
For  fixed  or  convecting  particles,  the  stability  of  the  particle  wakes  in  a  Ismirwr  bountlary  layer  was  analyzed  with  consideration  given  to  the  wake  position, 
width  and  the  velocity  deficit.  This  provides  the  initial  conditions  for  the  development  of  a  wave  packet  getKrated  by  a  particle  in  a  boundary  layer.  The 
developnwnt  of  the  wave  packets  generated  by  the  particles  is  then  performed.  This  report  describes  the  first  items  in  the  analysis  and  the  remainder  are 
presented  in  a  doctoral  thesis. 
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Abstract 


An  experimental  and  analytical  study  of  aspects  of  transition  induced  by  disturbances  from 
spherical  particles  in  laminar  boundary  layers  is  discussed.  The  generation  of  turbulent  wedges 
by  fixed  spherical  particles  in  a  laminar  boundary  layer  on  or  near  the  surface  of  a  flat  plate  is 
considered  experimentally  using  flow  visualization  with  fluorescent  dye  and  laser  Doppler 
velocimetry.  Turbulent  spots  generated  by  freely  convecting  spherical  particles  that  are  released 
in  the  freestream  to  fall  into  a  flat  plate  laminar  boundary  layer  and  impact  the  plate  are  also 
discussed.  A  combination  of  dye  flow  visualization  and  a  video  based  particle  tracking  technique 
was  used  to  study  the  convecting  particle  problem.  Although  the  Reynolds  number  at  the  critical 
condition  for  turbulent  wedge  generation  by  fixed  particles  and  turbulent  spot  generation  by 
convecting  particles  are  similar,  transition  in  these  two  situations  appears  to  be  fundamentally 
different.  The  development  of  a  turbulent  wedge  near  the  critical  condition  is  a  relatively  gradual 
process.  In  contrast,  turbulent  spots  form  relatively  quickly  after  the  convecting  particles  enter 
the  boundary  layer  and  impact  the  plate.  Turbulent  wedge  formation  downstream  of  a  fixed 
particle  results  ft-om  the  destabilization  of  the  near  wall  flow  by  the  vortical  structures  shed  into 
particle  wake.  This  shedding  process  is  dominated  by  periodically  shed  loop  shaped  hairpin 
vortices.  Observation  of  subharmonic  oscillations  at  and  *4  of  Uiis  shedding  frequency  suggest 
that  a  chaotic  route  to  turbulence  by  a  series  of  period  doubling  bifurcations  is  possible. 
Convecting  particle  wakes  have  little,  if  any,  affect  on  turbulent  spot  generation.  The 
disturbances  to  the  flow  near  the  wall  due  to  the  particle  impact  rapidly  develop  and  spots  begin 
to  emerge  at  the  particle  in  all  cases. 

The  analysis  first  predicts  the  trajectories  of  the  convecting  particles.  The  development  of  wave 
packets  from  disturbances  at  different  locations  is  considered.  For  fixed  or  convecting  particles, 
the  stability  of  the  particle  wakes  in  a  laminar  boundary  layer  was  analyzed  with  consideration 
given  to  the  wake  position,  width  and  the  velocity  deficit.  This  provides  the  initial  conditions 
for  the  development  of  a  wave  packet  generated  by  a  particle  in  a  boundary  layer.  The 
development  of  the  wave  packets  generated  by  the  particles  is  then  performed.  This  report 
describes  the  first  items  in  the  analysis  and  the  remainder  are  presented  in  a  doctoral  thesis. 
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Chapter  1 


Transition  Induced  by  a  Fixed  Spherical  Particle 
in  a  Laminar  Boundary  Layer 


Abstract:  The  development  of  the  wake  of  a  fixed  spherical  particle  mounted  on  the  surface  or 
elevated  above  the  surface  in  a  laminar  boundary  layer  and  the  resultant  downstream  formation 
of  a  turbulent  wedge  has  been  studied  experimentally.  Hairpin  vortices  are  shed  periodically  into 
the  wakes  of  these  particles  in  all  critical  cases  and  a  horseshoe  vortex  is  established  when  the 
particle  is  sufficiently  close  to  the  surface.  Elevated  particle  wake  data  indicate  the  horseshoe 
vortex  does  not  have  a  significant  effect  on  turbulent  wedge  formation.  The  legs  of  the  hairpin 
vortices  play  a  significant  role  in  the  breakdown  of  the  flow.  Periodicity  develops  in  these  wakes 
at  'A  and  %  of  the  hairpin  vortex  shedding  frequency  and  a  chaotic  route  to  turbulence  by  ? 
cascade  of  period  doubling  bifurcations  is  possible. 


1.1  Introduction 

Boundary  layer  transition  from  laminar  to  turbulent  flow  is  an  important  aspect  of  fluid 
dynamics  in  numerous  engineering  applications.  Natural  transition  in  laminar  boundary  layers 
involves  the  selective  amplification  of  small  disturbances  leading  to  the  development  of  two 
dimensional  instability  waves.  Secondary  instabilities  and  three  dimensionality  develops  and  the 
flow  rapidly  breaks  down  and  becomes  turbulent.  Alternate  routes  to  turbulence  involving  bypass 
transition  are  possible.  Bypass  transition  includes  transition  processes  that  bypass  the  2D  stage 
of  natural  transition  and,  typically,  involves  disturbances  that  are  initially  larger  and  more  three 
dimensional  than  in  the  case  of  natural  transition.  Laminar  flow  control  is  an  attractive 
technology  for  reducing  skin  friction  and  flow  noise.  Laminar  flow  can  be  maintained  to 
unnaturally  large  Reynolds  numbers  by  a  variety  of  methods  such  as  surface  suction,  surface 
heatingOn  liquids),  and  body  shaping.  The  degradation  of  laminar  flow  control  performance  by 
sources  of  disturbance  in  the  operating  environment  often  involves  bypass  transition  and  defeats 
this  technology. 

Lauchle,  et  al.,(1986)  evaluated  how  sized  particles  seeded  into  the  freestream  degraded 
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the  performance  of  a  laminar  flow  control  heated  body  in  a  water  tunnel.  Transition  Reynolds 
numbers  as  large  as  37  x  10^  were  achieved  in  a  clean  flow.  However,  freestream  particles 
decreased  the  transition  Reynolds  number  substantially.  No  minimum  particle  diameter  threshold 
for  performance  degradation  was  observed  but  larger  particles  had  a  more  substantial  effect.  In 
a  similar  heated  body  laminar  flow  control  experiment  examining  particle  degradation  of 
performance  by  Ladd  and  Hendricks(198S),  the  estimated  critical  particle  diameter  Reynolds 
numbers  based  on  the  velocity  in  the  undisturbed  flow  one  particle  diameter  above  the  surface 
at  the  location  turbulent  spots  were  observed  corresponds  i^proximately  to  values  reported  of  the 
critical  Reynolds  number  by  Hall(1967)  tor  fixed  spherical  particles  in  laminar  pipe  flow.  It  was 
speculated  by  Ladd  and  Hendricks(198S)  that  the  particles  impact  the  body  and  momentarily 
stick. 

Ha]l(1967)  used  flow  visualization  to  establish  conditions  resulting  in  turbulent  wedge 
formation  downstream  of  a  fixed  spherical  particle  mounted  on  or  near  the  wall  in  a  developing 
laminar  pipe  flow.  When  a  particle  was  positioned  near  the  wall,  the  effective  diameter  was  the 
distance  from  the  wall  to  the  top  of  the  particle,  k.  The  critical  Reynolds  number,  Re^  varied 
from  585  to  665  for  these  near  surface  cases.  The  particle  diameter  Reynolds  number  is  based 
on  the  velocity  at  the  top  of  the  particle  in  the  undisturbed  flow  throughout  this  ch^ter  of  the 
report.  Hali(1967)  cited  the  data  of  Magarvey  and  Bishop(1961)  to  suppon  the  hypothesis  that 
the  shear  layer  shed  from  the  these  fixed  spheres  is  turbulent  at  the  critical  condition  and  this 
turbulence  seeds  a  critical  region  of  the  boundary  layer  below  y/6  =  0.4  to  cause  the  formation 
of  a  turbulent  wedge.  Magarvey  and  Bishop(l%l)  observed  transition  occurs  at  Re^  »  600  in 
the  wake  of  a  sphere  in  a  uniform  freestream.  However,  these  results  differ  with  other  work 
cited  in  a  survey  of  sphere  wake  data  by  Sakamoto  and  Haniu(1990). 

Sakamoto  and  Haniu(1990)  describe  the  intricate  evolution  of  the  wake  of  a  sphere  in  a 
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uniform  flow  with  increasing  Reynolds  number.  Periodic  shedding  of  hairpin  vortices  begins  at 
Re^  =  300  -  400.  The  heads  of  these  hairpin  vortices  are  open  round  loops  and  not  tight 
hairpins.  This  is  the  case  for  all  of  the  vortices  described  as  hairpins  below.  The  vortex 
shedding  orientation  is  random  but  remains  fixed  at  these  low  Re^.  Irregularity  of  the  shedding 
frequency  and  side  to  side  motion  of  the  hairpin  vortices  occurs  intermittently  for  420  <  Re^  < 
480.  As  Re^  is  increased  above  Re^  =  480,  irregular  periodic  shedding  and  slow  irregular 
rotation  of  the  shedding  orientation  is  observed.  Above  Rej  =  650,  the  Strouhal  number  based 
on  the  hairpin  vortex  shedding  frequency  drops  from  0.21  to  0.19.  Above  Re^  =  800,  smaller 
vortex  loops  appear  in  the  shear  layer  shed  from  the  sphere  and  the  shedding  frequency  Strouhal 
number  bifurcates  into  an  upper  and  a  lower  frequency  mode.  Transition  from  a  laminar  to  a 
turbulent  wake  is  said  to  occur  with  the  bifurcation  at  Re^  =  800. 

Mochizuki(l961)  studied  the  flowfield  created  by  spherical  particles  on  the  surface  in  a 
laminar  boundary  layer  with  flow  visualization.  The  structures  in  the  wake  were  observed  over 
a  range  of  Reynolds  numbers.  As  with  Hall(1967),  a  turbulent  wedge  formed  downstream  of  the 
particle  when  the  Re^  reached  a  critical  value  and  the  wedge  moved  to  the  sphere  as  the  Reynolds 
number  was  increased.  The  critical  Reynolds  increased  from  approximately  Re^  ^  =  670  to 
1000  as  the  Reynolds  number  based  on  the  distance  of  the  particle  from  the  leading  edge.  Re,, 
and  freestream  velocity  was  decreased  from  220,000  to  below  50,000.  Based  on  these  data  and 
Hairs(1967)  data,  the  dependence  of  Re^  on  Re,  is  small  above  Re,  =  100,000.  It  was 
speculated  that  the  hairpin  vortex  heads  do  not  affect  turbulent  wedge  formation  and  that  the 
horseshoe  vortex  formed  by  the  flow  separation  upstream  of  the  sphere  and  the  legs  of  the  hairpin 
vortices  play  a  key  role  in  the  transition  process. 

Acarlar  and  Smith(1987)  and  Klebanoff,  Cleveland,  and  Tidstrom(1992)  examined  the 
wake  of  a  hemispherical  bump  in  a  flat  plate  laminar  boundary  layer.  The  interaction  between 
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the  horseshoe  and  hairpin  vortices  and  between  successive  hairpin  vortices  was  believed  lead  to 
the  development  of  secondary  vortices  and,  subsequently,  to  turbulent  wedge  formation. 

1.2  Objective 

Quantitative  information  in  the  wakes  of  surface  mounted  spheres  and  hemispheres  that 
generate  a  turbulent  wedge  is  limited  and  the  mechanisms  leading  to  development  of  a  turbulent 
wedge  are  poorly  understood.  The  case  of  a  fixed  elevated  sphere  in  a  laminar  boundary  layer 
has  hardly  been  considered.  The  objective  of  the  present  fixed  spherical  particle  study  was  to 
investigate  the  evolution  of  the  wakes  of  fixed  spherical  particles  on  or  near  the  surface  in  a 
laminar  boundary  layer  for  subcritical  to  supercritical  conditions  to  gain  a  greater  understanding 
of  the  details  of  these  wake  flows  and  the  processes  leading  to  formation  of  a  turbulent  wedge. 
Another  objective  of  these  fixed  particle  turbulent  wedge  data  deals  comparison  to  the  convecting 
particle  results  in  Chapter  2  of  this  report.  It  was  desired  to  determine  whether  these  phenomena 
of  turbulent  wedge  and  spot  formation  by  spherical  particles  appears  to  be  fundamentally  the  same 
or  whether  they  are  distinctly  different.  This  last  point  is  discussed  in  Chapter  2. 

1.3  Experimental  Apparatus  and  Approach 

Experiments  were  conducted  in  a  low  speed  closed  loop  open  channel.  The  laminar 
boundary  layer  developed  on  a  flat  plate  mounted  above  the  floor  of  this  channel  was  used  for 
this  transition  study.  Figure  1.1  is  diagram  of  the  channel  and  of  the  test  plate.  The  plenum  of 
the  channel  is  supplied  through  a  diffuser  containing  series  of  7  perforated  plates  which  aid  in 
breaking  up  large  structures.  This  is  followed  by  a  IS  cm  thick  section  of  3.2  mm  cell  size 
honeycomb  and  the  1:5.88  contraction  nozzle.  The  flat  plate  leading  edge  was  a  12:1  ellipse. 
The  channel  test  section  dimensions  are  a  75.9  cm  span  and  a  3.66  m  length  with  a  50  cm  depth. 
The  plate  span  is  2cm  less  than  the  channel  span  with  a  1  cm  gap  at  the  sidewalls.  This  sidewall 
gap  could  be  adjusted  with  movable  slats  at  the  bottom  of  the  plate  to  affect  a  bleed  with  which 
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sidewall  contamination  effects  might  be  controlled.  These  slats  were  set  so  that  the  gaps  were 
closed  for  the  present  study.  Observations  were  limited  to  the  upstream  portion  of  the  test  section 
to  avoid  sidewall  contamination.  The  sidewalls  of  the  test  section  were  glass  for  optical  access. 

Flow  visualization  was  achieved  by  seeping  a  concentrated  solution  of  fluorescein 
disodium  salt,  a  fluorescent  dye,  into  the  flow.  Wallcer(1987)  provides  details  about  the 
performance  of  the  fluorescein.  An  inclined  slit  in  the  surface  of  the  plate  located  upstream  of 
the  fixed  particle  location  and  a  strategically  placed  wand  located  in  the  plenum  were  used  to  add 
the  dye  to  the  flow.  Two  ultraviolet  ’black’  lights  emitted  enough  radiation  in  the  range  of  the 
absorption  spectrum  of  the  dye,  which  peaks  at  490  nm  wavelength,  to  produce  fluorescence 
effectively  without  lighting  up  the  background  or  the  plate  surface.  The  flow  visualizations  were 
videotaped  for  later  analysis. 

Velocity  profile  measurements  were  made  with  a  one  component  Helium-Neon  laser 
based  fringe  mode  forward  scatter  LDV  system  with  frequency  shift  in  one  of  the  two  laser 
beams.  A  TSI  model  1980b  digital  frequency  counter  and  a  PC  computer  were  used  to  process 
the  Doppler  signals  and  sample  LDV  data.  Vincent(1993)  provides  further  details  about  the 
instrumentation  and  the  facility. 

Velocity  profiles  of  the  streamwise  component  were  measured  on  the  centerplane 
upstream  of  a  surface  mounted  particle  in  the  horseshoe  vortex  system  that  forms  when  the 
approach  flow  separates.  The  wakes  of  surface  mounted  and  elevated  particles  were  surveyed 
from  X  =  1  to  as  far  as  X  =  240  radii  downstream  of  the  particle.  X  is  the  streamwise  distance 
from  the  particle  normalized  with  the  particle  radius.  Y  is  the  vertical  distance  from  the  plate 
similarly  normalized  with  the  particle  radius.  Near  critical  cases  generating  a  turbulent  wedge 
were  emphasized  and  representative  cases  are  discussed  herein.  Multiple  sets  of  LDV  data  were 
taken  for  estimating  velocity  fluctuation  power  spectra  using  the  Discretized  Lag  Product  method 
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discussed  by  Mayo,  et  al.(1974).  Surveys  in  the  wakes  of  spheres  in  the  uniform  freestream  were 
made  for  comparison  purposes.  Vincent(1993)  presents  additional  velocity  survey  dau  not 
included  in  this  report.  This  includes  limited  vertical  component  profiles  and  Reynolds  shear 
stress  profiles. 

Spherical  particles  were  mounted  in  position  on  or  above  the  surface  at  1  =  0.838  m 
downstream  of  the  leading  edge  of  the  test  plate  using  0.2S  mm  diameter  wire.  The  wire 
extended  from  the  plate  to  the  particle  when  the  particles  were  on  or  near  the  surface.  This  was 
the  case  for  the  two  cases  detailed  in  this  Ch^ter  and  the  wire  alone  was  not  shedding  a  vonex 
street  due  to  the  low  local  Reynolds  number.  The  wire  extended  through  the  particle  to  a  support 
above  the  free  surface  for  particles  further  from  the  surface.  A  Questar  QMl  long  range 
microscope  with  a  graduated  reticle  eyepiece  and  a  translation  stage  with  a  0.025  mm  position 
readout  was  used  to  evaluate  the  position  of  the  particles  relative  to  the  wall. 

1.4  Results 

Figure  1.2  compares  the  measured  mean  velocity  profiles  on  the  centerline  of  the  channel 
to  the  Blasius  profile  for  the  freestream  velocities  used  of  the  cases  detailed  in  the  discussion 
below.  The  99%  thickness  is  denoted  as  5  and  y  is  the  distance  from  the  surface.  There  is  a 
slight  favorable  pressure  gradient  in  the  channel  due  to  an  inclined  flap  at  the  trailing  edge  of  the 
test  plate  and  a  section  of  honeycomb  above  the  working  surface  at  the  exit  of  the  test  section. 
The  tail  flap  and  honeycomb  created  a  slight  negative  angle  of  attack  at  the  nose  of  the  test  plate 
to  avoid  unsteady  separation  on  the  working  surface.  The  freestream  turbulence  intensity 
measured  by  the  LDV  system  was  0.39%.  The  actual  intensity  is  probably  less  since  the  He-Ne 
based  LDV  signals  processed  by  a  digital  counter  are  expected  to  be  noisy  relative  to,  say,  a  hot¬ 
wire. 
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1.4.1  Flow  Visualization  with  a  Surface  Mounted  Sphere 

A  particle  sufficiently  near  the  surface  separates  the  approach  flow  and  a  horseshoe  vortex 
system  is  formed.  The  horseshoe  vortex  wraps  around  the  particle  and  the  legs  of  the  vortex 
form  two  streamwise  oriented  counter-rotating  vortices  that  extend  downstream  into  the  wake  of 
the  particle.  When  dye  is  seeped  into  the  flow  at  the  wall  upstream  of  the  particle,  the  separation 
line  at  the  edge  of  the  horseshoe  vortex  is  well  defined  because  none  of  the  surface  injected  dye 
is  entrained  past  this  separation  line  into  the  horseshoe  vortex.  The  dye  laden  surface  flow  is  not 
lifted  but  is  diverted  spanwise  around  the  horseshoe  vortex  system.  A  distinct  keyhole  shaped 
dye  shadow  results  and  defines  the  perimeter  of  the  horseshoe  vortex  well  downstream  into  the 
wake.  The  dye  on  the  surface  is  drawn  towards  the  centerplane  of  the  particle  wake  slightly 
downstream  of  the  particle  as  the  flow  converges  in  the  region  of  hairpin  vortex  formation.  Dye 
seeped  into  the  flow  from  a  wand  in  the  plenum  is  needed  to  visualize  the  horseshoe  and  hairpin 
vortices  in  the  wake  of  the  particles. 

The  present  estimates  of  the  critical  Reynolds  number  for  surface  mounted  particles  are 
consistent  with  the  values  and  trends  observed  by  Mochizuki(1961)  and  Hall(1967).  A  continuous 
turbulent  wedge  was  observed  in  the  present  study  at  Re;,crit  =  630  to  700  with  Re,  =  129,000. 
Rea  ^  is  based  on  the  particle  diameter  and  the  velocity  at  the  top  of  the  particle  in  the 
undisturbed  flow.  Mochizuki(1961)  observed  a  wedge  forming  at  Re^  cnt  =  68S  and  Re,  = 
122,000.  Defining  when  a  continuous  turbulent  wedge  forms  is  subject  to  some  personal  bias 
because  the  wake  is  irregular  appearing  before  the  wedge  is  clearly  present.  Below  Re^  =  480 
to  SOO,  the  hairpin  vortices  do  not  lift. 

Frames  of  the  video  recorded  flow  visualizations  were  digitized  with  a  frame  grabber  and 
contrast  enhanced  because  the  images  are  dim.  Figure  1.3  is  a  series  of  these  images  showing 
side  views  at  successive  streamwise  positions  for  a  d  =  6.33  mm  diameter  sphere  on  the  surface 
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with  U,  =  14.3  cm/s,  d/5  =  0.55,  Re,  =  129,000  and  Re^  ==  810.  A  turbulent  wedge  develops 
downstream  of  the  sphere  for  this  case.  The  displacement  thickness  Reynolds  number  at  the 
sphere  position  in  the  undisturbed  flow  was  Re^.  -  540. 

A  slight  side  to  side  oscillation  of  the  hairpin  vortex  heads  develops  in  the  wake.  This 
cannot  be  detected  in  the  side  views  of  figure  1.3  but  is  observed  from  above  the  plate.  Each 
successive  vortex  head  moves  away  from  the  vertical  centerplane  in  a  direction  opposite  to  the 
adjacent  heads.  A  similar  motion  is  seen  in  the  wakes  of  spheres  in  a  uniform  flow,  see 
Sakamoto  and  Haniu(1990). 

Figure  1.3a  shows  the  region  from  tqiproximately  X  =  0  to  80.  The  flow  moves  from 
left  to  right.  The  hairpin  vortices  begin  forming  near  X  =  5  and  begins  to  rise  near  X  =  10. 
Hie  legs  of  the  hairpin  vortices  descend  to  near  the  surface  at  X  =  25.  The  vortex  legs  draw 
closely  together  as  they  extend  to  the  surface.  Depending  on  the  position  of  filaments  of  dye, 
narrow  fountains  or  jets  of  dye  can  occasionally  be  seen  lifting  vertically  between  the  legs 
immediately  upstream  of  the  vortex  heads.  The  vortex  heads  are  round  open  loops.  Table  1.1 
lists  the  height  of  the  hairpin  tops  at  various  streamwise  locations  and  brief  observations  from  the 
flow  visualizations.  By  X  =  50,  the  tops  of  the  heads  have  lifted  past  the  edge  of  the  boundary 
layer  in  the  undisturbed  flow. 

The  legs  of  the  hairpin  vortices  away  from  the  wall  are  inclined  at  steeper  angles  as  they 
move  downstream.  The  inclination  of  the  straight  segment  of  the  legs  increases  from  30°  to  35° 
initially  after  formation  to  45°  in  figure  1 .3a.  Fluctuations  can  be  detected  near  the  wall  in  figure 
1.3a  after  the  vortex  legs  descend  to  the  wall.  This  development  continues  in  figure  1.3b,  X  = 
64  to  130.  The  fluctuations  near  the  surface  where  the  vortex  legs  are  present  are  more  turbulent 
appearing  in  the  close-up  view  further  downstream  in  figure  1.3c,  X  =  176  to  205.  Given  that 
the  dye  filaments  retain  a  history  of  fluid  motions,  the  flow  may  not  be  as  irregular  as  it  appears. 


However,  the  velocity  profile  data  presented  below  support  that  an  increased  level  of  fluctuation 
is  developing  near  the  wall.  The  hairpin  heads  appear  to  separate  from  the  legs  of  the  vortices 
and  begin  to  look  like  isolated  vorticity  rings.  The  flow  becomes  increasingly  irregular  at  the 
wall  and  the  heads  of  the  hairpins  appear  to  have  little  to  do  with  this  irregularity  as  speculated 
previously  by  Mochizuki(1961)  and  Acarlar  and  Smith(1987). 

Figure  1.3d  shows  that  the  third  vortex  head  from  the  right  side  is  descending  towards 
the  surface  and  a  bulge  of  turbulent  appearing  fluid  is  apparent  immediately  upstream  of  this 
downward  moving  vortex.  A  definite  periodic  pattern  develops  with  every  fourth  hairpin  drawn 
towards  the  surface. 

1.4.2  Flow  Visualization  with  an  Elevated  Sphere 

Moving  the  sphere  away  from  the  surface  results  in  a  wake  flow  that  develops  with  the 
same  general  features  as  the  surface  mounted  case  but  there  are  differences.  The  wall  has  a 
stabilizing  influence  on  the  hairpin  vortices;  the  vortices  are  shed  with  fixed  orientation  so  that 
the  heads  lift  directly  away  from  the  plate  and  not  with  a  random  irregular  orientation.  For  the 
critical  conditions  of  interest,  the  hairpins  vortices  move  noticeably  from  side  to  side  with  some 
slight  rotation  about  an  axis  through  the  sphere.  This  rotation  increases  and  becomes  irregular 
as  the  distance  between  the  sphere  and  the  surface  is  increased.  The  sphere  in  a  uniform  flow 
condition  is  approached  with  further  elevation  of  the  particle.  Also,  the  horseshoe  vortex  that 
forms  upstream  of  the  surface  mounted  particle  diminishes  in  size  and  strength  as  the  sphere  is 
moved  away  from  the  wall  until  it  is  no  longer  evident  in  the  surface  flow  visualizations.  Exactly 
when  the  horseshoe  vortex  system  is  lost  as  the  particle  is  elevated  was  not  systematically  studied 
but  for  the  cases  considered  the  vortex  diminishes  quickly  as  the  gap  between  the  surface  and  the 
bottom  of  the  sphere  was  increased  beyond  O.lSd  and  the  horseshoe  vortex  is  not  detected  with 
the  flow  visualization  with  gaps  larger  than  0.2Sd  to  0.3Sd. 
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Figure  1.4  shows  side  views  of  the  wake  of  a  d  =  4.72  mm  diameter  sphere  located 
above  the  plate  at  the  same  streamwise  position  and  freestream  velocity  as  the  surface  mounted 
case  in  figure  1.3.  Re*  »  690  for  this  elevated  particle  and  a  continuous  turbulent  wedge  forms 
near  X  =  100  in  the  wake.  The  bottom  of  the  sphere  is  at  y/6  =  0.33  and  the  top  is  at  k/6  = 
0.74  with  d/5  =<  0.41;  k  is  the  distance  from  the  wall  to  the  top  of  the  particle.  The  gap  between 
the  bottom  of  the  sphere  and  the  plate  was  0.8d  and  the  horseshoe  vortex  could  not  be  detected. 

Figure  1.4a  shows  that  the  hairpin  vortex  heads  lift  away  from  the  surface  more  quickly 
than  with  the  surface  mounted  case.  The  legs  of  the  hairpin  vortices  extend  to  the  near  the 
surface  at  X  =  30  which  is  further  downstream  than  the  surface  mounted  case.  The  hairpin 
heads  grow  in  size  and  the  inclination  of  the  legs  varies  further  downstream  in  figure  1.4b.  As 
with  the  surface  mounted  case,  the  legs  of  the  hairpin  vortices  are  inclined  at  30°  to  35°  initially 
but  this  angle  does  not  increase  quickly  to  45°.  Well  away  from  the  wall,  the  legs  kink  and  are 
nearly  vertical.  The  heads  of  the  vortices  are  further  away  from  the  wall  than  with  the  surface 
mounted  case  considered  above  and  this  draws  the  legs  up  into  a  region  of  low  mean  shear. 
Table  1.2  summarizes  observations  from  the  flow  visualization. 

Fluctuations  near  the  wall  begin  to  move  away  from  the  surface  and  some  of  the  hairpin 
heads  move  down  towards  the  wall  with  a  repeating  pattern;  every  fourth  vortex  head  in  figure 
1.4c  is  slightly  lower  than  the  other  vortices.  Further  downstream,  the  lower  frequency 
periodicity  observed  with  the  surface  mounted  case  at  '4  the  shedding  frequency  is  readily 
apparent  for  this  elevated  particle  case.  The  more  obvious  appearance  of  this  14  shedding 
frequency  cycle  compared  to  the  surface  mounted  case  may  result  from  the  diminution  of  the 
stabilizing  influence  of  the  wall  on  wake  as  the  particle  is  elevated  above  wall. 

Figures  1.4c  and  1.4d  show  that  every  fourth  hairpin  deforms  and  moves  towards  the 
wall  and  turbulent  bulges  of  lifting  fluid  tqipear  immediately  upstream  of  the  descending  vortex 
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heads.  Details  such  as  the  shape,  size,  and  intensity  distribution  of  the  dye  in  the  vortex  heads 
repeats  with  amazing  regularity  in  the  region  covered  by  figure  1.4d.  The  hairpin  heads  blend 
together  between  the  lifting  turbulent  bulges  of  fluid  and  it  appears  that  the  outward  motion  of 
the  bulges  is  accompanied  by  an  inward  rush  of  fluid  that  pulls  the  vortex  heads  together  and 
towards  the  surface.  The  cycle  is  similar  to  the  cycle  of  bursting  and  sweeps  observed  in  the 
near  wall  region  of  a  turbulent  boundary  layer  by  Corino  and  Brodkey(1969)  but  these  turbulent 
bulges  are  also  reminiscent  of  the  low  Reynolds  numbers  turbulent  boundary  layer  flow 
visualizations  by  Falco(1977)  showing  similar  coherent  motions  in  the  outer  flow. 

The  transition  of  the  wake  of  cylinders  was  studied  by  Kamiadakis  and 
Triantafyllou(1992)  and  Tomboulides,  et  al.(1992),  with  a  Navier-Stokes  equation  simulation. 
Period  doubling  bifurcations  were  observed  in  these  simulations  after  the  vortex  street  develops 
secondary  instabilities.  The  flow  alternates  between  two  preferred  states  and  the  limit  cycle  in 
phase  space  branches  into  two  connected  cycles.  The  first  period  doubling  is  the  result  of  half 
wavelength  shifts  of  the  spanwise  modulation  of  the  flow.  A  cascade  of  period  doubling 
bifurcations  with  increased  Re^  occurs  and  leads,  ^patently,  to  a  chaotic  route  to  turbulence. 

Although  the  wakes  of  spheres  cannot  support  a  spanwise  flow  modulation  as  with  a 
cylinder,  the  orientation  of  the  hairpin  vortex  head  about  the  axis  of  the  sphere  could  be 
considered  as  an  analogous  component  or  degree  of  freedom.  The  side  to  side  rotational  motion 
of  the  hairpin  vortex  heads  at  half  the  shedding  frequency  observed  in  the  present  flow 
visualizations  is  a  case  of  the  flow  alternating  between  two  preferred  states.  The  regular  motion 
of  every  forth  hairpin  vortex  head  towards  the  surface  could  also  be  another  period  doubling 
bifurcation.  A  r^id  breakdown  of  the  flow  follows  as  suggested  by  the  emergence  of  the 
turbulent  bulges. 
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1.4J  Horseshoe  Vortex  Velocity  Profiles 

Velocity  profiles  were  measured  at  10  streamwise  locations  upstream  of  a  surface 
mounted  sphere  from  X  =  -l.Ol  to  -3.56  for  Re^  =  700  with  a  particle  diameter,  d  =  3.91,  a 
freestream  velocity,  U,  =  23.3  cm/s,  and  Re^,  =  710  where  6'  is  the  displacement  thickness. 
Hie  flow  visualizations  in  figure  1.2  used  a  larger  particle  and  a  lower  freestream  velocity.  The 
conditions  used  for  these  velocity  profile  measurements  result  in  d/d  =  0.46  which  nearly 
matches  the  elevated  particle  case  in  figure  1.4.  The  Re,,  values  are  nearly  matched  too. 
However,  Re^.  is  lower  for  the  elevated  case,  Rej.  =  540. 

The  intensity  of  the  fluctuations  in  the  horseshoe  vortex  region  upstream  of  the  sphere, 
based  on  the  freestream  velocity,  reaches  a  peak  value  of  9.2%  at  y/k  =  0.5  or  Y  =  1 
immediately  upstream  of  the  sphere  at  X  =  -1.01  as  shown  in  figure  1.5  where  k  =  d  for  this 
surface  mounted  case.  This  measurement  location  is  close  to  the  stagnation  point.  This  intensity 
peak  near  y/k  =  0.5  decreases  to  under  3%  at  X  =  -1.26.  These  fluctuations  probably  result 
from  the  periodic  shedding  and  not  a  horseshoe  vortex  oscillation.  The  dye  flow  visualizations 
show  that  the  incoming  flow  is  divided  and  drawn  into  the  horseshoe  vortex  for  y/k  <  0.5  or 
passes  around  the  sphere  for  y/k  >  0.5  similar  to  the  findings  of  Mochizuk?^.  Therefore,  the 
largest  expected  horseshoe  vortex  diameter  is  comparable  to  the  radius  of  the  sphere. 

The  velocity  fluctuation  intensity  normalized  with  U.  increases  near  the  surface  in  the 
profile  at  X  =  -1.01  in  figure  1.5.  This  near  wall  peak  in  fluctuation  intensity  increases  moving 
upstream  away  from  the  sphere  with  a  maximum  values  from  5%  to  5.5%  at  X  =  -1 .26  to  -1 .51 . 
By  X  =  -2.53,  the  fluctuation  intensities  near  the  surface  are  similar  to  the  levels  near  the  wall 
with  no  particle.  Mean  reverse  flow  near  the  surface  was  detected  as  far  upstream  as  X  =  -2.02 
with  a  maximum  observed  magnitude  of  0.06U.  at  X  =  -1.51.  This  is  the  general  location  of 
maximum  fluctuation  intensity  near  the  wall. 
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1.4.4  Mean  Velocity  Profiie  Evolution  in  the  Wakes 

The  development  of  the  wake  and  turbulent  wedge  generated  by  the  surface  mounted 
sphere  (d  =  3.91  mm,  U,  =  23.3  cm/s,  Re^  =  700,  d/6  =  0.46)  was  evaluated  by  measuring 
velocity  profiles  on  the  centerplane  at  X  =  1.01,  3,  5,  10,  20,  50,  100,  150,  200,  and  240. 
Select  mean  profiles  are  shown  in  figure  1.6.  Inflections  due  to  the  presence  of  the  sphere  and 
hairpin  vortices  decay  and  move  away  from  the  surface  as  the  profile  gradually  takes  a  form 
similar  to  a  turbulent  boundary  layer.  Figure  1.7  is  a  similar  set  of  mean  profiles  for  the  elevated 
particle  case  (d=4.72  mm,  U,  =  14.3  cm/s,  Re^  =*  690,  y/8  =  0.41,  k/8  =  .74)  that  was 
considered  in  figure  1 .4.  The  dashed  lines  on  the  figures  indicates  the  location  of  the  top  and 
the  bottom  of  the  sphere.  Figures  1.6  and  1.7  show  data  that  are  averaged  over  many  vortex 
shedding  periods  without  regard  to  the  phase  of  the  shedding  cycle.  As  a  result,  these  mean 
profiles  and  the  profiles  discussed  below  may  differ  substantially  from  instantaneous  profiles. 

The  influence  of  the  hairpin  vortices  and  the  momentum  deficit  of  the  sphere  on  the 
velocity  profiles  is  more  pronounced  for  the  elevated  case  in  figure  1 .7  than  the  surface  mounted 
case  in  figure  1.6.  The  wake  of  the  elevated  sphere  initially  has  a  local  velocity  minima  near  the 
particle  center  as  is  the  case  for  a  sphere  in  the  freestream.  The  hairpin  vortices  extend  further 
from  the  surface  for  the  elevated  case  than  the  surface  mounted  case  and  this  influences  the 
velocity  profiles  further  away  from  the  wall  and  for  a  greater  streamwise  distance  than  the  wall 
mounted  case.  The  elevated  sphere  profile  at  X  =  100  is  the  solid  line  on  figure  1.7  and  exhibits 
a  step  near  U/U.  =  0.7.  Apparently,  during  the  point  by  point  sequence  of  data  taking  during 
the  survey,  the  turbulent  wedge  moved  from  upstream  of  this  location  to  downstream  of  it.  As 
a  result,  the  profile  is  turbulent  like  below  U/U.  =  0.7  and  laminar  like  at  the  larger  velocity 
values  further  from  the  wall. 

Figure  1.8  shows  the  development  of  the  shape  factor,  H  =  6*/6,  for  the  two  cases 
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considered  in  figures  1.6  and  1.7  and  a  subcritical  surface  mounted  particle  case  that  shed 
hairpins  th' .  do  not  lift.  The  Blasius  value  is  H  =  2.59.  The  shape  factor  trend  is  similar  to  that 
observed  by  Schubauer  and  Klebanoff(1955)  for  a  transitional  boundary  layer.  The  elevated 
sphere  wake  develops  more  quickly  than  the  surface  mounted  case  even  though  Re^  and  d/d  are 
similar  for  the  two  cases  and  Rej.  is  larger  for  the  surface  mounted  case. 

An  estimate  of  the  friction  velocity,  u*,  was  made  at  the  most  downstream  stations  using 
a  procedure  similar  to  Chang  and  Blackwelder(1990).  Spalding’s  law  of  the  wall  is  fit  to  the 
profile  data  using  u*  as  the  free  parameter.  Only  data  near  the  wall  are  used  to  avoid  the  hairpin 
distortions  to  the  mean  profile  away  from  the  wall.  The  resultant  mean  profiles  are  plotted  in 
inner  variaoles  where  u*  =  U/u*  and  y*  =  yu*/v.  These  profiles  are  compared  with 
Spalding’s(1961)  iaw  of  the  wall  in  figure  1.9.  The  profiles  have  a  log  region  and  are  turbulent- 
like  in  most  respects.  The  wake  region  of  the  elevated  particle  boundary  layer  profile  is  still 
appreciably  distorted.  The  surface  mounted  particle  velocity  profile  at  X  =  240  is  replotted  in 
outer  variables  and  compared  with  a  low  Reynolds  number  turbulent  profile  from  Purteli,  et 
al.(1981)  in  figure  1.10.  The  data  below  y/5*  =  3  suggests  the  flow  is  still  developing. 
Considering  the  periodic  structures  observed  in  the  flow  visualizations  in  figures  1.3  and  1.4  at 
these  downstream  stations,  the  similarity  with  turbulent  profiles  is  noteworthy. 

1.4.5  Velocity  Fluctuation  Intensity  Profiles 

Figures  1.11  and  1.12  show  the  evolution  of  ffeestream  velocity  normalized  turbulence 
intensity  profiles  at  select  locations  for  the  surface  mounted  and  elevated  spheres  considered  in 
figures  1 .6  and  1 .7.  There  are  two  peaks  in  the  fluctuation  intensity  profiles  beginning  at  X  = 
10  and  continuing  downstream  until  the  outer  peak  is  lost.  The  outer  peak  corresponds  to  the 
heads  of  the  hairpin  vortices  and  moves  away  from  the  surface  and  decays  with  increasing 
streamwise  distance.  This  is  similar  to  results  for  a  hemispherical  bump  wake  in  Acarlar  and 
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Sinith(1987)  and  KJebanoff,  et  al.(l992).  For  the  surface  mounted  case,  the  inner  and  outer 
intensity  peaks  are  about  the  same  at  X  =  50  but  this  matching  of  inner  and  outer  intensity  peaks 
is  estimated  to  occur  near  X  =  25  to  30  for  the  elevated  particle  case.  At  X  =  100,  both  cases 
still  have  an  outer  turbulence  intensity  peak  even  though  the  mean  profiles  show  no  obvious 
feature  at  X  =  100.  The  elevated  particle  intensity  peak  is  almost  twice  as  far  from  the  wall  as 
the  surface  mounted  case.  The  outer  peak  is  not  distinct  at  X  =  200  in  either  case  and  is  beyond 
the  upper  limit  of  the  survey  for  the  elevated  particle. 

For  the  surface  mounted  case,  the  outer  turbulence  intensity  peak  is  30%  to  32%  at  X 
=  I  and  3  and  it  decreases  to  28%  and  18  %  at  X  =  5  and  10.  For  the  elevated  particle  case, 
the  outer  peak  intensities  near  the  particle  are  slightly  larger;  47%,  37%,  24.4%,  and  15.6%  at 
X  =  1,  3,  5,  and  10,  respectively.  Although  these  intensity  levels  are  high,  they  meet  the  10% 
to  18%  levels  mentioned  by  Morkovin(1969)  for  self  sustaining  turbulence,  and  the  mean  profiles 
undergo  inflections  associated  the  hairpin  heads,  the  distance  of  these  disturbances  from  the 
surface,  their  frequency  and  their  range  of  scales  are  such  that  instability  mechanisms  are  not 
directly  activated  as  suggested  by  their  decay. 

For  comparison  purposes,  the  wakes  of  spheres  in  a  freestream  were  surveyed  at  various 
Reynolds  numbers  including  Re,,  =  700,  760,  and  1000  at  X  =  1,  2,  3,  5,  10  and  20.  The 
conditions  Re,,  =  700,  760,  and  1000  span  a  range  comparable  to  the  near  critical  cases 
considered  and,  therefore,  are  of  some  interest.  Two  distinct  intensity  peaks  due  to  the  shear 
layer  shed  from  the  sphere  were  observed  off  the  centerline  of  these  freestream  spheres  at  Re^ 
=  700  and  above.  The  intensity  levels  were  35%  to  45%  at  X  =  1,2,  and  3  in  this  Re  range. 
At  X  =  5,  the  peaks  move  towards  the  centerline  and  decay  to  levels  near  25%.  By  X  =  10, 
the  intensity  maxima  associated  with  the  hairpin  vortices  are  flat  across  the  wake  at  levels  near 
12%  and  this  decreases  to  near  6%  at  X  =  20.  This  flat  intensity  profile  does  not  occur  for  the 
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near  surface  spheres  because  the  orientation  of  the  shed  vortices  is  fixed.  The  intensities 
observed  with  the  surface  mounted  particle  are  initially  lower  than  observed  with  the  freestream 
particles  but  the  decay  in  the  intensities  near  and  just  downstream  of  where  the  hairpins  form, 
X  =  10,  is  less  than  with  the  freestream  particles.  The  intensity  peaks  for  the  elevated  particle 
associated  with  the  hairpin  vortex  heads  mimic  those  seen  with  the  freestream  particles.  In  a 
limited  sense,  the  presence  of  the  wall  and  boundary  layer  have  little  apparent  effect  on  the 
hairpin  vortex  heads  other  than  to  suppress  the  randomness  of  shedding  orientation. 

The  near  wall  peak  in  the  fluctuation  intensity  first  appears  at  X  =  10  for  both  the 
surface  mounted  and  the  elevated  particle  cases.  This  intensity  peak  slowly  grows  in  vertical 
extent  and  level  with  increasing  downstream  distance  as  shown  in  figures  1.11  and  1.12.  The 
intensities  gradually  approach  turbulent  boundary  layer  levels  but  the  peak  does  not  move  to  the 
wall  until  downstream  of  X  =  100  in  both  cases.  This  is  where  the  turbulent  wedge  develops. 
At  X  =  10,  the  near  wall  peak  is  8%  and  is  located  at  y/d  =  0.27  for  the  surface  mounted  case. 
The  level  at  X  =  10  for  the  elevated  particle  case  was  5%  and  the  peak  is  located  at  y/k  =  0.25, 
approximately.  Acarlar  and  Smith(1987)  observe  a  similar  double  peaked  intensity  profile  as  near 
as  X  =  4  downstream  of  a  hemispherical  bump  but  the  near  wall  peak  is  further  from  the  wall 
relative  to  the  distance  to  the  top  of  the  particle,  at  Y  =  0.5. 

Activity  near  the  wall  upstream  of  where  the  near  wall  turbulence  intensity  maxima  are 
first  detected,  X  =  10,  does  not  appear  to  be  a  significant  tac;.  At  X  =  1,3,  and  5  the 
fluctuation  intensity  profiles  decrease  monotonically  as  the  wall  i'-  approached  for  both  surface 
mounted  and  elevated  cases  and  no  features  of  obvious  importance  are  observed.  The  fluctuation 
levels  near  the  wall  at  these  locations  are  different  for  the  two  cases  because  the  elevated  particle 
data  near  the  wall  were  taken  below  the  particle  and  outside  of  the  particle  wake.  However,  both 
cases  still  develop  a  fluctuation  peak  near  the  wall  by  X  =  10  that  is  not  seen  at  X  =  5.  The 
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flow  visualizations  indicate  X  =  10  is  about  where  the  hairpin  heads  have  clearly  formed  and 
begin  lifting  but  the  legs  of  the  vortices  still  extend  back  to  near  the  top  of  the  spheres  rather  than 
towards  the  wall.  As  noted  for  a  hemispherical  bump,  the  hairpins  form  downstream  of  a  region 
of  flow  convergence  behind  the  particles.  The  circulation  induced  by  the  vortex  tubes  forming 
the  hairpin  lifts  low  speed  fluid  near  the  wall  and  can  create  periodically  varying  inflectional 
velocity  profiles.  This  lifted  low  speed  fluid  and  the  pressure  gradients  in  the  hairpin  formation 
region  may  be  key  to  the  emergence  of  the  disturbances  at  X  =  10. 

The  movement  of  the  particle  away  from  the  wall  weakens  the  horseshoe  vortex  and  the 
flow  visualizations  indicate  that  the  horseshoe  vortex  was  not  present  for  the  elevated  case 
considered.  Given  that  the  development  of  the  surface  mounted  and  elevated  particle  wakes  near 
the  wall  is  generally  the  same,  the  role  of  the  horseshoe  vortex  in  the  wake  development  and 
turbulent  wedge  formation  process  is  probably  not  significant. 

1.4.6  Hairpin  Vortex  Shedding  Frequencies 

Velocity  fluctuation  power  spectra  were  determined  for  surface  mounted  spheres  on  a 
limited  basis.  The  LDV  probe  volume  was  positioned  near  hairpin  vortex  head  passage  locations 
in  an  effort  to  determine  the  hairpin  vortex  shedding  frequency.  The  data  were  taken  at  Y  =  3 
and  4  and  at  X  =  20  and  36  with  3.21  mm  and  3.91  mm  diameter  particles  and  between  X  = 
10  and  16  at  Y  =  3  with  a  S.Sl  mm  diameter  particle  over  a  range  of  Reynolds  numbers.  At 
the  lower  Reynolds  numbers  consider,  Re^  <  800,  the  spectra  were  flat  versus  frequency,  except 
for  narrow  spikes  due  to  the  shedding  frequency.  At  X  =  20  and  36  and  Rq,  =  1000  to  1500, 
the  spectra  began  to  broaden  and  fill  near  the  shedding  frequencies.  This  suggests  the  shedding 
process  is  becoming  more  complex  or  that  variations  in  the  interactions  between  the  wake 
structures  occur  to  cause  cycle  to  cycle  phase  shifts. 

Figure  1.13  plots  Strouhal  numbers,  St  =  fd/U,,  of  distinct  peaks  in  the  spectra  of 
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surface  mounted  spheres  versus  Re^.  These  data  are  clustered  into  linear  bands  illustrated  by  the 
line  segments  on  figure  1.13.  The  two  bands  at  lower  Strouhal  values  are  similar  to  the  high  and 
low  frequency  shedding  modes  that  appear  with  a  sphere  in  a  uniform  flow  described  in  Sakamoto 
and  Haniu(1990).  The  high  frequency  mode  results  are  at  slightly  greater  St  than  for  a 
hemispherical  bump  where  St  is  based  on  the  radius  of  the  bump  in  Acarlar  and  Smith(1987). 
The  two  higher  frequency  bands  on  figure  1.13  are,  apparently,  harmonics  of  the  high  frequency 
mode.  There  are  only  two  data  points  in  the  highest  frequency  band  so  it  is  difficult  to  determine 
much  about  these  data. 

The  14  shedding  frequency  subharmonic  development  of  turbulent  bulges,  observed  well 
downstream  for  the  near  critical  Reynolds  number  cases  considered  above,  are  not  detected  in 
figure  1.13  even  at  higher  Re^.  The  low  Strouhal  number  band  of  data  is  at  or  near  the 
shedding  frequency  and  may  be  a  subharmonic  associated  with  the  side  to  side  motion  of  the 
vortex  heads.  The  low  frequency  shedding  mode  Strouhal  number  data  in  Sakamoto  and 
Haniu(1990)  for  a  sphere  in  a  uniform  flow  deviate  from  this  type  of  behavior  and  are  nearly 
constant  with  values  close  to  0.2  in  range  600  <  Re^  <  2000. 

1.5  Summary  and  Conclusions 

The  development  of  the  wakes  of  surface  mounted  and  elevated  particles  in  a  laminar 
boundary  layer  was  studied  using  flow  visualization  and  laser  Doppler  velocimetry.  The  wake 
of  a  surface  mounted  particle  that  generated  a  turbulent  wedge  developed  more  gradually  than  an 
elevated  particle  wake  at  matched  conditions.  However,  the  elevated  particle  wake  development 
was  generally  similar  to  the  surface  mounted  particle  wake.  This  suggests  that  horseshoe  vortex 
that  forms  with  the  surface  mounted  particle  does  not  play  a  significant  role  in  the  wake 
development  and  turbulent  wedge  formation  process.  A  thin  layer  of  irregular  fluctuating  fluid 
emerges  close  to  the  wall  near  X  =  10  for  both  surface  mounted  and  elevated  particle  cases. 
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This  X  position  is  where  the  hairpin  vortices  have  formed  and  begin  lifting  fluid  away  from  the 
surface.  This  fluctuating  layer  near  the  wall  grows  in  the  vertical  extent  of  the  flow  it  includes 
and  its  intensity  level  grows  with  increasing  streamwise  distance.  Eventually  the  intensity  levels 
approach  turbulent  boundary  layer  levels  and  the  fluctuation  intensity  peak  moves  to  the  wall  as 
the  turbulent  wedge  develops.  The  legs  of  the  hairpin  vortices  appear  to  play  a  role  in  the 
development  of  secondary  structures  in  this  layer. 

The  hairpin  vortices  shed  by  the  spheres  oscillate  side  to  side  at  one  half  of  the  shedding 
frequency  and  periodic  ejections  of  turbulent  bulges  occur  further  downstream  at  one  quarter  of 
the  shedding  frequency.  These  subharmonic  motions  were  more  pronounced  with  the  elevated 
particle  case.  This  may  indicate  a  relaxation  of  the  effects  of  the  wall  on  the  particle  wake  as  the 
particle  is  elevated  above  the  surface  and  may  explain  the  more  rapid  development  of  the  elevated 
particle  wake.  The  development  of  these  periodicities  at  double  and  quadruple  the  shedding 
period  suggests  that  a  chaotic  route  to  turbulence  by  a  cascade  of  period  doubling  bifurcations 
is  possible. 
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Figure  l.I  Diagram  of  open  channel  with  test  plate  installed. 


Fig.  1.2  Mean  velocity  proxies  in  undisturbed  boundary  layer  compared  with  the  Blasius 
profile. 
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Fig.  1.3  (a)  top  and  Cb)  bottom.  See  next  page  for  caption. 


Fig.  1.4  (a)  top  and  (b)  bottom.  See  next  page  for  caption. 


Fig.  1.4  Side  views  of  elevated  panicle  wake.  Re^  =  690.  k/5  =  0.74,  and  d/5  =  0.41. 
(c)  top  and  (d)  bottom,  (a)  .X  =  0  -  65.  (b)  X  =  107  -  !50.  (c)  X  =  172  -  215.  (d)  X 
=  306  -  339. 
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Fig.  1.6  Mean  velocity  profiles  in  the  surface  mounted  sphere  wake. 


o 

X=  10 

X 

X  =  20 

X=  100 

A 

X  =  200 

0.0  0.2  0.4  0.6  0.8 

U/U- 

Fig.  1.7  Mean  velocity  profiles  in  the  wake  of  the  elevated  sphere. 
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Fig.  1.8  Shape  factor,  H  =  h’l6,  for  a  subcritical  surface  mounted  sphere  and 
the  surface  mounted  and  elevated  cases  of  interest. 
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Fig.  1.9  Comparison  of  mean  profiles  in  inner  variables  with  Spalding ’s(  1961) 
law  of  the  wall. 
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Fig.  1.10  Comparison  of  the  surface  mounted  sphere  wake  mean  velocity  profile 
at  X  =  240  with  turbulent  boundary  layer  data  from  Purtell,  et  al(1981). 
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Fig.  1.11  Turbulence  intensity  profiles  in  the  wake  of  the  surface  mounted 
sphere. 
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Fig.  1.12  Turbulence  intensity  profiles  in  the  wake  of  the  elevated  sphere. 
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Fig.  1.13  Hairpin  vortex  shedding  frequency  Strouhal  numbers  versus  Reynolds  number. 
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Table  1.1  Observations  based  on  the  surface  mounted  sphere  flow  visualizations. 


X 

Top  of 
Hairpin,  y/5 

Comments 

1 

— 

no  hairpin  vortices 

3 

— 

no  hairpin  vortices 

5 

— 

vortex  roll  up  starting 

10 

0.46  -  0.50 

hairpin  vortices  form 

20 

0.80 

hairpin  vortices  rising,  boundary  layer 
orderly  near  wall 

50 

1.16 

disturbances  near  the  wall,  hairpin  legs 
reattach  to  surface  and  flatten 

100 

1.40 

hairpin  legs  beginning  to  break  up 

150 

1.60 

intermittent  turbulent  activity  near  wall 

240 

1.60 

turbulent  bulges,  every  4th  hairpin  head 
moves  toward  plate 

Table  1.2  Observations  based  on  the  elevated  sphere  flow  visualizations. 


X 

Top  of 
Hairpin,  y/6 

Comments 

1 

— 

no  horseshoe  vortex,  no  hairpin  vortices 

3 

— 

1.5d  separation  bubble  attached  to  sphere 

5 

0.87 

hairpin  heads  begin  to  form  and  rise 

10 

1.13 

hairpin  heads  rising,  legs  at  height  of 
bottom  of  sphere  and  moving  toward  plate 

20 

1.30 

hairpin  legs  approach  plate,  heads  rising 

50 

1.74 

hairpin  legs  reattach  to  plate  surface  at  X  = 
30,  hairpin  heads  rising 

100 

2.17 

boundary  layer  irregular  near  wall, 
turbulent  bursts  evident,  every  4th  vortex 
head  is  low 

200 

2.18 

boundary  layer  chaotic,  turbulence  bulges 
affect  hairpin  heads,  periodic  pattern 
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Chapter  2 


Turbulent  Spot  Formation  by  Freely  Convecting 
Spherical  Particles  in  a  Laminar  Boundary  Layer 


Abstract:  The  generation  of  turbulents  spots  from  disturbances  created  by  spherical  particles  that  fall 
freely  into  a  laminar  boundary  layer  was  studied  experimentally  using  flow  visualization  and  a  video 
based  particle  tracking  technique.  Turbulent  spots  formed  when  the  particle  velocity  Reynolds 
number  exceeded  700.  Although  turbulent  wedges  form  at  a  similar  critical  Reynolds  number  in  the 
wake  of  fixed  surface  mounted  spheres,  the  two  processes  are  not  the  same.  Turbulent  spots,  when 
they  occur,  emerge  soon  after  the  disturbance  generating  particle  impacts  the  surface  and  the  spots 
emerge  in  the  immediate  vincinity  of  the  particles.  Turbulent  wedges  form  downstream  of  fixed 
spherical  particles  after  a  gradual  development  and  growth  of  disturbances  near  the  wall  in  the  wake 
of  the  particles.  Features  such  as  periodically  shed  hairpin  vortices,  observed  with  fixed  particles,  are 
not  observed  in  the  freely  convecting  particle  cases. 


2.1  Introduction 

Laminar  flow  control  technology  has  the  potential  to  reduce  skin  friction  and  flow  noise  for 
marine  and  aeronautical  applications.  Laminar  flow  can  be  maintained  to  unnaturally  large  Reynolds 
numbers  by  a  variety  of  methods  such  as  surface  suction,  surface  heating(in  liquids),  and  body 
shaping  to  maintain  a  favorable  pressure  gradient.  However,  the  degradation  of  laminar  flow  control 
performance  on  a  body  by  the  generation  of  turbulent  spots  due  to  surface  roughness,  by  disturbances 
generated  by  particles  suspended  in  the  freestream  flow,  and  by  freestream  turbulence  defeats  this 
technology  in  anything  but  a  clean  fluid  dynamic  environment  typical  to  a  laboratory.  The 
practicality  of  this  technology  for  field  applications  is  limited  as  a  result. 

The  generation  of  the  turbulent  patches  on  a  laminar  flow  surface  from  environmental 
disturbances  often  involves  bypass  transition.  Bypass  transition  refers  to  processes  that  bypass  the 
linear  amplification  of  2D  waves  of  natural  transition  and,  as  discussed  by  Morkovin(l969),  may 
include  numerous  paths  to  turbulence.  Similar  bypass  transition  problems  can  affect  turbomachinery 
applications.  Mayle(1991)  discusses  how  bypass  transition  induced  by  a  variety  of  disturbances  can 
be  a  significant  factor  affecting  the  performance  and  design  of  turbomachinery. 
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Lauchle  and  Gurney(1984)  quantified  the  performance  of  a  laminar  flow  control  heated  body 
in  a  water  tunnel.  Lauchle,  et  al.(1986)  studied  how  seeding  the  freestream  with  sized  particles 
degraded  this  performance.  Transition  Reynolds  numbers  as  large  as  37  x  10*,  based  on  the  surface 
arc  length  on  the  body  and  the  approach  flow  velocity,  were  achieved  in  a  clean  flow.  However, 
freestream  particles  decreased  the  transition  Reynolds  number  substantially.  No  minimum  particle 
diameter  threshold  for  performance  degradation  was  observed  but  larger  particles  had  a  more 
substantial  effect.  Details  of  the  mechanisms  leading  to  breakdown  of  laminar  flow  were  not 
determined.  In  a  similar  heated  body  laminar  flow  experiment,  Ladd  and  Hendricks(1985)  observed 
that  the  Reynolds  numbers  based  on  the  diameter  of  particles  seeded  into  the  freestream  and  the 
estimated  velocity  a  distance  one  particle  diameter  above  the  surface  at  the  streamwise  position  where 
turbulent  spots  formed  on  the  body,  Rq^  erit>  corresponds  approximately  to  the  Reynolds  numbers  that 
Hall(1967)  observed  transition  induced  by  fixed  spherical  particles  on  or  near  the  surface  in  a 
developing  laminar  pipe  flow.  Ladd  and  Hendricks(198S)  speculated  that  the  particles  producing  the 
turbulent  spots  impact  the  body  and  momentarily  stick. 

Chen,  Goland  and  Reshotko(1979)  modelled  the  ciqiture  of  freestream  suspended  particles  into 
the  boundary  layer  of  a  laminar  flow  body  and  predicted  the  rate  of  turbulent  spot  generation.  Spot 
generation  rate  prediction  was  based  on  an  empirically  based  particle  size  criteria  since  the  actual 
mechanisms  are  not  known.  For  a  given  particle  size  distribution  and  a  fixed  specific  gravity,  the 
turbulent  spot  generation  rate  could  be  predicted  as  a  function  of  the  freestream  velocity. 

Vortical  structures  in  the  wake  of  fixed  particles  mounted  on  or  near  the  surface  in  a  laminar 
boundary  layer  developed  gradually  into  a  turbulent  wedge  at  the  critical  Reynolds  number  as 
discussed  in  Chapter  1.  Although  strong  disturbances  are  created  in  the  form  of  hairpin  vortices  in 
the  near  wake  of  these  fixed  particles,  these  disturbances  decay.  Tliese  observations  support  the 
speculation  by  Mochizuki(1961)  that  the  heads  of  the  hairpin  vortices  have  little  to  do  with  turbulent 
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wedge  formation  in  the  wakes  of  surface  mounted  spheres.  The  results  in  Chapter  1  support  that  the 
horseshoe  vortex  formed  by  the  flow  separation  upstream  of  surface  mounted  spheres  did  not  have  an 
effect  on  turbulent  wedge  formation.  The  legs  of  the  hairpin  vortices  and  the  disturbances  near  the 
wall  appear  to  cause  the  turbulent  wedge  to  develop.  Whether  these  fixed  particle  experiments  are 
relevant  to  the  case  of  a  freestream  particle  crossing  streamlines  and  impacting  the  surface  of  a 
laminar  flow  control  body  has  not  been  determined.  Ladd  and  ’/iendricks(198S)  suggestion  that  the 
particles  momentarily  stick  at  the  surface  to  cause  a  turbulent  spot  suggests  that  this  is  the  case. 

A  small  part  of  the  experiments  performed  by  Hall(1967)  in  a  developing  laminar  pipe  flow 
included  pulling  tethered  spherical  particles  into  the  wall.  Particles  that  were  too  small  to  generate  a 
turbulent  wedge  when  fixed  at  the  wall  could  cause  a  turbulent  spot  when  pulled  to  the  wall  with 
sufficient  velocity. 

2.2  Objectives 

A  series  of  experiments  were  performed  to  provide  a  qualitative  evaluation  of  the  types  of 
disturbances  generated  near  the  surface  by  freely  falling  particles  that  enter  a  laminar  boundary  layer 
and  impact  the  surface  and  to  quantitatively  determine  the  conditions  resulting  in  turbulent  spot 
formation.  The  present  effort  is  not  intended  to  be  an  exhaustive  and  definitive  study  but,  given  the 
lack  of  knowledge  of  the  phenomena  involved,  is  a  first  step  towards  gaining  a  detailed  understanding 
of  the  mechanisms  of  spot  generation  by  impacting  suspended  particles. 

2.3  Experimental  Apparatus  and  Procedures 

The  experiments  were  conducted  in  the  low  speed,  closed  loop,  open  channel  facility 
described  in  Chapter  1  and  in  Vincent(1993)  using  the  laminar  boundary  layer  that  develops  on  a  test 
plate  mounted  above  the  floor  of  the  channel  test  section.  Particles  were  released  from  rest  to  fall 
from  a  location  14  cm  above  the  plate  and  22. S  cm  downstream  of  the  leading  edge  of  the  plate.  At 
the  fixed  freestream  velocity  used  for  all  of  the  experiments,  U.  =  14.3  cm/s,  the  particles  impact 
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the  plate  at  locations  between  30  and  60  cm  downstream  of  the  leading  edge.  This  is  as  much  as 
twice  the  distance  the  fixed  particles  were  located  downstream  from  the  leading  edge  in  Lhapter  1 
but,  as  mentioned  in  Chapter  1,  the  effect  on  the  critical  Reynolad  number  of  increasing  Re,  in  the 
range  used  is  small. 

The  probe  holding  a  given  particle  at  the  release  point  was  a  2. 1  mm  outside  diameter 
stainless  steel  tube  with  a  beveled  tip  to  support  the  particle.  The  probe  was  clamped  in  position  and 
pierced  the  free  surface  of  the  open  channel  flow.  A  vacuum  was  drawn  on  a  250  ml  flask  with  a 
hand  pump  and  this  evacuated  flask  was  connected  to  the  support  probe  by  a  piece  of  flexible  plastic 
tubing.  A  valve  was  opened  to  relieve  the  slight  vacuum  holding  the  particle  on  the  end  of  the  probe 
to  allow  the  particle  to  fall  away  from  the  probe.  The  cylindrical  probe  Reynolds  number  was  323 
and  the  probe  sheds  a  vortex  street.  This  vortex  street  could  affect  the  particle  paths  after  release. 
However,  the  aim  of  the  experiment  is  to  achieve  a  range  of  conditions  at  impact  and  observe  the 
results.  Therefore,  the  potential  effects  the  probe  may  have  on  the  particle  path  are  not  necessarily 
significant  complications.  The  particle  paths  were  smooth  almost  linear  arcs  with  no  obvious  probe 
effects  in  most  cases. 

The  vertical  component  of  velocity  of  the  particles  at  impact  with  the  test  plate  is  a  function 
of  the  particle  diameter  and  the  particle  specific  gravity  for  the  fixed  release  location  and  freestream 
velocity.  Polystyrene,  nylon,  and  delrin  spheres  3.1  mm  to  7.9  mm  in  diameters  with  nominal 
specific  gravities  of  l.OS,  1.17,  and  1.41  were  used  to  achieve  a  range  of  impact  conditions.  The 
sphericity  and  diameter  tolerance  of  the  spheres  was  ±0.001”  and  ±0.002".  The  diameters  of 
individual  spheres  used  were  measured  with  a  micrometer. 

Visualization  of  the  disturbances  created  at  the  wall  during  and  after  particle  impact  was 
achieved  by  seeping  a  fluorescent  dye  solution,  fluorescein  disodium  salt,  through  a  narrow  inclined 
slot  at  the  surface  of  the  plate.  Walker(1987)  discusses  aspects  of  using  fluorescene.  Visualizations 
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were  done  for  a  range  of  conditions  and  videotaped  from  above  the  plate  to  determine  which  cases 
produced  a  turbulent  spot.  The  whole  process  was  repeated  and  videotaped  from  the  side  to 
determine  velocity  and  pathline  histories.  Since  all  of  the  particles  paths  are  influenced  by  their 
unsteady  wake,  some  variation  in  the  conditions  at  impact  between  the  top  and  side  view  replications 
is  likely  and  some  laditude  for  this  variation  must  be  granted  when  interpreting  the  results. 

The  sideviews  of  the  particle  paths  were  analyzed  frame  by  frame.  After  digitizing  the 
particle  positions  in  a  frame  by  frame  play  back  of  the  video  tape  to  obtain  the  paths,  velocity 
histories  were  obtained  by  differentiating  the  paths.  Since  the  process  of  differentiation  of  data  is 
suceptable  to  errors  from  uncertainty  in  the  pathline  data,  a  cubic  spline  least  squares  algorithm  with 
hing  was  applied  to  the  pathline  data  to  estimate  the  velocity  histories. 

2.  Results 

Table  2.1  lists  information  pertaining  to  the  particles  used,  the  impact  conditions  and  location, 
the  whether  or  not  a  turbulent  spot  is  generated  after  impact.  Estimates  of  the  terminal  vertical 
velocities  are  listed  for  each  particle.  These  terminal  velocities  are  estimated  using  an  empirical  drag 
coefficient  given  in  White(1974)  and  presented  below.  One  trial  for  each  case  and  viewing 
orientation  was  recorded.  Analysis  of  the  video  was  done  frame  by  frame  with  human  analysis  at 
eac  ^p.  This  process  was  slow  and  the  number  of  trials  was  limited  as  a  result. 

2.4.1  Flow  Visualization  Results 

In  addition  to  determining  if  a  turbulent  spot  is  formed,  the  flow  visualizations  show  details  of 
the  disturbances  at  the  wall  throughout  the  process  of  a  particle  impacting  and  then  moving  along  the 
plate.  Typical  illustrative  examples  are  presented  in  this  section.  In  all  cases,  the  particles  rolled  on 
the  surface  after  settling  from  impact  but  it  was  not  determined  whether  there  was  slipage  with  this 
rolling  motion.  The  majority  of  the  particles  considered  bounced  at  impact.  Figure  2. 1  shows  a 
sequence  of  images  taken  with  a  S.S6  mm  diameter  polystyrene  particle.  The  first  image,  figure 
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2.1a,  is  approximately  at  impact.  The  flow  moves  from  left  to  right  and  d/5  =  0.62  and  d/5*  ==  1.85 
at  impact.  This  case  corresponds  to  particle  #16  in  table  2.1.  The  Reynolds  number  based  on  the 
two  component  velocity  magnitude  of  the  particle,  Uu,,  at  impact  is  between  1000  and  1100. 

A  scar  in  the  dye  begins  forming  at  the  wall  just  prior  to  impact  and  near  impact,  figure  2.6a, 
dye  has  been  gathered  up  into  two  stream* vise  oriented  tubes  or  streaks  with  the  particle  in  between. 
The  source  of  light  is  from  the  bottom  of  the  figures  and  the  particle  can  be  detected  from  the  shadow 
it  leaves.  The  particle  convects  faster  after  impact  than  the  dye  laden  flow  near  the  surface  so  that  the 
dye  streaks  aligned  with  the  freestream  flow  are  left  on  the  upstream  side  of  the  particle  in  its  wake. 
After  impact,  the  U  component  particle  velocity  decreases  from  near  U«  to  near  0.75U.  in  time 
period  of  about  0.25s.  The  velocity  in  the  undisturbed  boundary  layer  is  approximately  0.85U,  at 
the  top  of  the  particle,  assuming  the  particle  is  on  the  surface.  The  resultant  slip  velocity  at  the  top  of 
this  moving  particle  is  only  O.IU.  and  the  resulting  slip  velocity  Reynolds  number  is  only  about  85. 
Spanwise  oscillations  in  the  dye  streaks  are  readily  apparent  on  the  upstream  side  of  the  particle  in 
figures  2.1a  near  impact  and  2.1b  0.66s  after  impact.  Initially,  these  oscillations  may  have  something 
to  do  with  the  one  noticeable  bounce  after  impact  but  the  oscillations  continue  after  impact  and  could 
indicate  a  disturbance  amplification  stage.  Although  the  velocity  of  the  particle  relative  to  the  fluid  is 
low  at  the  top  of  the  particle,  this  is  not  the  case  near  the  wall  where  the  fluid  velocities  are  lower. 
Disturbances  on  the  downstream  side  of  the  particle  are  detected  in  figure  2.1b  and  these  develop 
quickly  into  what  begins  to  look  like  an  emerging  turbulent  spot  next  to  the  particle  in  figure  2.1c 
which  is  1.33  s  after  impact.  The  delayed  appearance  of  what  looks  like  instabilities  followed  by  a 
rapid  breakdown  of  the  flow  does  fit  the  general  scenario  described  by  Breuer  and  Landahl(1990)  for 
localized  strong  disturbances.  The  spot  is  fully  formed  2.0  s  after  impact  and  figure  2.  Id  shows  this 
spot.  Notice  the  streaks  left  behind  the  spot  on  the  upstream  side  of  the  particle  straighten  out  in 
figures  2.1c  and  2.  Id  but  spanwise  oscillation  of  the  streaks  are  visible  in  the  turbulent  spot.  This  is 
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can  be  seen  in  previous  flow  visualizations  of  turbulent  spots,  see  Elder(l960)  for  example.  In  this 
particular  case,  the  particle  convects  with  the  spot  at  the  trailing  edge  of  the  spot. 

Smaller  particles  were  observed  to  lag  the  spot  they  generated  and  did  not  generate  additional 
spots  or  a  turbulent  wedge.  Larger  particles  moved  with  the  spots,  often  staying  near  the  leading 
edge  of  the  spots.  The  difference  in  this  convective  velocity,  presumably,  is  that  the  particles  convect 
with  a  velocity  that  increases  with  size  due  to  the  effect  of  the  velocity  gradient. 

Figures  2.2a,  b,  and  c  show  the  flow  visualization  of  a  subcritical  3.97  mm  diameter 
polystyrene  particle  at  impact  and  1  s  and  2  s  after  impact.  The  particle  leaves  a  dark  dye  shadow  on 
its  upstream  side  and  two  bright  dye  streaks  extend  downstream  of  the  particle.  This  case 
corresponds  to  sphere  #18  on  table  2.1.  At  impact,  d/5  =  0.42  and  d/6*  =  1.26  and  the  velocity  at 
the  top  of  the  particle  in  the  undisturbed  boundary  layer  is  about  0.67U,.  The  U  velocity  of  the 
particle  slows  to,  roughly,  0.3U.  shortly  after  impact.  This  corresponds  to  a  velocity  below  y/6  = 
0.2  in  the  undisturbed  boundary  layer.  Since  the  dye  streaks  are  moving  faster  than  the  particle,  they 
are  probably  lifted  farther  from  the  wall  than  y/6  =  0.2.  The  slip  velocity  at  the  top  of  the  particle  is 
actually  larger  than  the  case  considered  in  figure  2.1  and  the  Reynolds  number  based  on  this  slip 
velocity  is  about  225.  At  this  Re^,  based  on  the  results  in  Chapter  1,  the  particle  is  probably  not  be 
shedding  hairpin  vortices  periodically  but  if  vortices  are  shed,  they  do  not  lift. 

Figure  2.3a  shows  a  d  =  6. 35  mm  polystyrene  particle,  #15  on  table  2.1,  just  after  estimated 
impact.  The  scar  formed  in  dye  layer  at  the  surface  is  wishbone  shaped  with  streaks  at  oblique 
angles.  Figure  2.8b  is  0.66  s  later  and  after  the  particle  has  landed  from  the  bounce  on  impact.  The 
half  angle  of  the  oblique  streaks  in  the  dye  is  25°.  This  angle  is  larger  than  the  half  angle  enclosing 
the  area  passed  over  by  a  growing  turbulent  spot.  The  typical  maximum  half  angle  subtended  by  a 
spot  from  its  virtual  origin  is  10°,  see  Wygnanski,  et  al.(1976).  The  dye  panerns  and  disturbances 
generated  by  the  double  impact  of  this  bouncing  larger  particle  are  complex  enough  that  more  subtle 
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features  that  might  be  present  such  as  the  spanwise  oscillation  of  the  streamwise  streaks  noticed  in  the 
case  for  figure  2.1  are  not  detected.  A  spot  begins  to  emerge  about  1  s  after  impact  and  figure  2.3c 
shows  the  emerging  spot  1.33  s  after  impact.  The  spot  might  be  considered  slightly  more  developed 
than  the  case  considered  in  figure  2.1.  However,  some  of  the  disturbances  in  figure  2.3c  at  the 
spanwise  edges  of  the  emerging  spot  appear  to  be  artifacts  of  the  impact  disturbance  and  not  part  of 
the  spot.  The  edges  of  what  looks  like  the  emerging  spot  are  within  a  10°  wedge  extending  from 
near  the  point  of  impact. 

The  time  required  after  impact  for  the  spot  to  begin  emerging  in  the  cases  considered  is  tU./6 
=«  IS  to  20  or  tU./6*  =:  45  to  60.  In  the  simulation  of  the  evolution  of  strong  localized  disturbances 
by  Breuer  and  Landahl(1990),  the  magnitude  of  u’/U.  was  0.05  after  tU./S*  =  43  but  rapid  growth 
was  not  observed  until  tU./5*  =117  and  after.  The  present  disturbances  are  not  what  was  simulated 
and  the  development  into  a  spot  is  rapid  in  comparison. 

The  dye  streaks  formed  at  various  stages  of  the  events  in  figures  2.1  -  2.3  could  be  the  result 
of  shear  layers,  waves,  or  tubes  of  vorticity  above  the  dye  layer  at  the  wall.  This  exact  nature  of 
these  disturbances  is  difficult  to  judge  because  the  dye  at  the  wall  does  not  show  structural  features 
that  may  be  present  away  from  the  wall  unless  the  dye  is  first  lifted  from  the  wall.  The  dye  is  not 
lifted  significantly  until  a  spot  begins  to  emerge.  The  blocking  effect  of  the  wall  on  the  fluid  that  the 
approaching  particle  displaces  may  result  in  the  initial  scar  formation  for  the  larger  faster  falling 
particles  such  as  shown  in  figAire  2.3a.  It  is  significant  that  the  dye  is  not  lifted  appreciably  until  the 
spot  forms.  Anything  as  vigorous  as  the  hairpin  vortices  associated  with  the  fixed  particle  case  in 
Chapter  1  is  not  likely.  The  trailing  wake  left  by  these  particles  as  they  fall  is  not  expected  to  have  a 
large  effect  near  the  wall  either.  The  polystyrene  particles  used  in  the  cases  considered  in  figures  2.1 
-2.3  were  all  in  near  equilibrium  with  the  U  component  of  velocity  for  much  of  the  freefall  time 
period  prior  to  impact.  The  trailing  wake  should  extend  above  the  particle  along  a  vertical  line  at 
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impact  and  convect  with  the  freestream  velocity.  The  spots  emerge  from  near  the  surface 
immediately  adjacent  to  the  particle  for  critical  cases. 

2.4.2  Impact  Conditions 

The  observed  trends  in  the  conditions  at  impact  are  generally  consistent  but  there  are 
variations.  These  variations  are  probably  due  to  the  unsteady  and  random  character  of  these  falling 
particle  wakes  with  increasing  particle  Reynolds  number.  The  wake  of  the  probe  near  the  release 
point  could  possibly  affect  the  particle  trajectory  and  the  resulting  conditions  at  impact.  Errors  in  the 
measured  values  are  possible  contributors  to  this  variation  of  the  impact  velocity  data  also.  The 
uncertainty  in  the  measured  particle  velocities  is  estimated  as  no  more  than  10%  of  the  freestream 
velocity.  Vincent(1993)  provides  an  error  analysis. 

To  estimate  the  scatter  typical  to  the  measurements  of  the  impact  velocities,  consider  the 
polystyrene  particle  data  in  table  2.1.  The  polystyrene  particles  come  into  near  equilibrium  with  the 
streamwise  velocity  of  the  freestream  flow  before  impacting  the  plate.  The  polystyrene  particles  are 
lighter  and  fall  more  slowly  than  the  other  materials.  Therefore,  they  accelerate  from  rest  to  the 
freestream  velocity  more  quickly  than  the  other  materials  and  they  have  more  time  to  approach 
equilibrium  before  impacting  the  plate  because  they  fall  more  slowly  than  the  denser  particles.  The 
six  impact  U  component  velocities  for  the  polystyrene  particles  are  all  near  the  freestream  velocity 
with  a  standard  deviation  from  the  freestream  velocity  of  4%  and  the  largest  deviation  is  7.7%.  Note 
that  the  seventh  polystyrene  particle,  #19  on  table  2.1  did  not  impact  the  plate  in  the  field  of  view. 

The  streamwise  velocity  at  impact  with  nylon  and  delrin  particles  shows  the  largest  variations 
in  results  with  some  values  low  relative  to  the  rest;  particles  1,  6,  7,  and  9  are  examples  of  this  type 
of  deviation.  In  all  cases,  the  vertical  velocity  at  impact  listed  on  table  2.1  is  less  than  the  terminal 
velocity  estimate  for  a  sphere  in  a  still  fluid.  The  smaller  lighter  particles  are  closest  to  terminal 
vertical  velocity. 
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2.4  J  Conditions  Generating  Turbulent  Spots 

A  velocity  associated  with  the  particle  at  impact  should  scale  with  the  size  of  the  resultant 
disturbance  and,  therefore,  correlate  with  turbulent  spot  generation.  Certainly,  if  the  particle  U 
velocity  were  in  near  equilibrium  with  the  freestream  flow  prior  to  impact,  the  vertical  velocity  at 
impact.  Via,,  might  be  expected  to  scale  well  with  the  resulting  disturbance  since  it  is  the  only  velocity 
component  out  of  equilibrium  with  the  flow.  Figure  2.4  plots  the  ratio  versus  Re^vim,  the 

Reynolds  number  based  on  the  vertical  impact  velocity.  The  solid  symbols  on  the  figure  did  not 
generate  a  turbulent  spot.  A  few  points  on  figure  2.4  do  not  fit  the  trends.  For  example,  particle  #17 
on  table  2.1  has  low  Vi„  and  Re^vim  values.  The  trend  in  the  data  indicates  that  the  critical  values  of 
RCdvim  increase  with  increasing  Vi„,/U,.  This  increase  in  the  critical  Reynolds  number  with 
disturbance  strength  was  not  expected.  Based  on  other  transitional  flow  situations,  it  is  expected  that 
a  threshold  critical  Reynolds  number  may  exist  as  with  boundary  layers  or  sphere  wakes  or  that  the 
critical  Reynolds  number  might  decrease  with  increasing  disturbance  strength  similar  to  roughness 
element  effects.  This  unexpected  behavior  suggests  that  Vi„  alone  does  not  scale  with  the  strength  of 
the  disturbances  leading  to  spot  formation.  It  is  possible  that  pathline  differences  affect  the  outcome 
but  this  view  is  not  favored.  The  heavier  particles  produce  the  larger  Vi„/U,  and  have  a  steeper 
impact  angle  than  the  smaller  lighter  particles.  The  difference  in  the  inclination  angle  for  the 
subcritical  conditions  nearest  to  the  critical  state  varies  from  27°  to  50°  going  from  polystyrene  to 
delrin. 

The  two  component  slip  velocity  based  on  the  difference  between  the  impact  velocity  and  the 
freestream  velocity  was  considered  to  correlate  critical  events  with  impact  conditions.  The  slip 
velocity  magnitude  is  AU  =  [(Ui^  -  U.)^  +  However,  the  U  component  of  the  slip  velocity 

is  not  large  for  the  near  critical  cases  of  interes;  and  the  resulting  trends  are  similar  to  those  in  figure 
2.4  with  Vi.. 
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The  Reynolds  numbers  based  on  the  slip  velocities,  Vi„,  or  AU,  characterize  Jie  state  of  the 
particle  wake  at  impact.  The  findings  of  fixed  spherical  particle  study  in  Chapter  I  suggest  that  even 
large  disturbances  generated  away  from  the  wall  but  in  the  boundary  layer  have  little,  if  any,  direct 
effect  on  the  generation  of  a  turbulent  wedge.  Disturbances  near  the  wall,  detected  by  fluctuation 
intensity  peaks  below  y/5  =  0.2S  in  Chapter  1,  appear  to  lead  to  flow  break  down  in  the  fixed 
particle  wakes.  This  is  supported  by  the  present  flow  visualizations  that  show  the  spots  emerge  from 
the  surface  and  this  appears  to  be  the  result  of  distubances  at  or  below  y/5  =  0.2.  There  is  no  clear 
reason  to  expect  that  the  wakes  trailing  a  convecting  particle  would  affect  turbulent  spot  generation. 

The  slip  velocities  mentioned  above  do  not  quantify  the  size  of  the  disturbance  the  particle 
generates  near  the  wall  at  impact.  The  combined  two  component  impact  velocity  magnitude,  = 
(Ui,B^  +  where  Ui,,  is  the  streamwise  component  of  the  particle  velocity  at  impact  and  not  a 

slip  velocity,  should  scale  well  with  the  size  of  the  disturbances  generated  near  the  wall  by  the 
particles  at  impact.  The  reasoning  is  that  near  wall  fluid  velocities  are  low  and  the  data  indicate  the 
particle  does  not  slow  noticeably  just  prior  to  impact.  Figure  2.5  shows  the  impact  velocity 
magnitude  versus  the  Reynolds  number  based  on  The  conditions  generating  turbulent  spots  near 
the  critical  condition  tend  to  cluster  but  this  is  difficult  to  judge  given  the  coarseness  of  the  data  grid 
achieved  with  the  various  diameter  particles  and  specific  gravities  that  were  used.  Actual  values  of 
Re^  are  listed  in  table  2.1.  The  sparse  data  matrix  and  the  uncertainty  in  the  values  means  that 
identifying  critical  Reynolds  number  trends  accurately  is  not  possible  unless  there  are  reasonably  large 
and  consistent  variations  over  the  range  of  conditions  and  this  is  not  the  case.  The  critical  Reynolds 
number  is  constant  within  the  ability  to  ascertain  the  trend.  Allowing  for  a  10%  uncertainty,  the 
rang<^  for  the  critical  condition  that  overlaps  for  the  three  materials  is  Re„  =  700  to  800.  It  appears 
that  particle  size  relative  to  the  boundary  layer  thickness  and  the  size  of  the  disturbance,  has 
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little  effect  and  Re^  is  constant.  This  is  supported  by  the  heated  body  experiment  by  Lauchle,  et 
al.(1986). 

The  influence  of  the  size  of  particles  seeded  into  the  ffeestream  on  the  performance  of  a 
laminar  flow  control  heated  body  was  considerd  in  experiments  by  Lauchle,  et  al(1986).  The 
Reynolds  number  based  on  the  approach  flow  velocity  and  particle  diameter  was  plotted  versus  d/6' 
where  5*  is  the  displacment  thickness  at  the  upstream  most  location  that  spots  emerged  on  the  body. 
These  critical  condition  data  collapsed  and  the  trend  is  shown  in  figure  2.6  which  is  reproduced  from 
Lauchle,  et  al.(1986).  Data  include  3  heating  conditions  from  27  to  75  kW  input  to  heat  the  body. 

In  the  present  experiment,  the  Reynolds  number  based  on  is  felt  to  be  equivalent  ot  the  approach 
flow  based  value.  The  d/5*  values  in  the  present  experiment  span  a  range  from  1.1  to  3.5, 
approximately,  with  the  conditions  straddling  the  critical  state  spaiming  d/5*  =  1.1  to  1.85.  In 
Lauchle,  et  al,(I986)  the  values  spanned  a  range  from  0.03  <  d/5*  <  0.7  and  it  was  observed  that 
for  d/5*  >  0.5,  the  critical  freestream  based  Reynolds  number  was  constant,  within  the  data  scatter, 
with  values  in  the  range  690  <  Revert  <  770,  see  figure  2.6.  As  noted  above,  the  estimated  critical 
Reynolds  number  in  the  current  experiment  is  similar. 

t 

The  observed  transition  Reynolds  number  for  these  convecting  particles  is  similar  to  the  Rea  ^ 
values  based  on  the  velocity  in  the  undisturbed  flow  at  the  top  of  the  particle  that  were  observed  with 
fixed  particles  in  Chapter  1  which  are  630  <  Rea  <  700.  The  transition  of  the  wake  of  spheres  in 
a  uniform  flow  occurs  at  Rea  near  800,  see  Sakamoto  and  Haniu(1990).  The  near  agreement  between 
these  three  sets  of  critical  Reynolds  numbers  may  be  fortuitous.  The  fixed  particle  transition 
discussed  in  Chapter  1  appears  to  depend  on  the  lifting  of  low  speed  fluid  by  the  vortices  in  the  wake 
and  the  interaction  of  the  legs  of  hairpin  vortices.  At  the  critical  condition  a  turbulent  wedge 
develops  slowly  as  the  disturbances  caused  by  the  fixed  particle  grow.  Transition  in  the  case  of  a 
sphere  in  a  uniform  ffeestream  involves  a  free  shear  layer  instability.  The  impacting  particles 
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generate  a  turbulent  spot  by  a  mechanism  that  does  not  appear  to  involve  either  hairpin  vortices  or  the 
instabilities  in  shear  layer  shed  into  the  wake  of  the  falling  particle. 

2.4.4  Pathline  and  Velocity  History  Data 

Although  the  pathline  and  velocity  history  data  prior  to  impact  may  not  directly  affect 
turbulent  spot  generation,  these  data  are  interesting  and  are  considered  below.  The  pathlines  traced 
by  the  particles  as  they  fell  to  the  surface  were  usually  mild  arcs  that  were  nearly  straight  lines.  The 
arcs  are  concave  downward  indicating  the  V  velocity  is  increasing  relative  to  the  U  velocity  as  is 
expected  for  a  particle  in  freefall.  Figure  2.7  contains  an  example  of  a  typical  particle  pathline, 
particle  ifS  on  table  2.1  which  is  a  d  =  3.94  mm  diameter  nylon  particle.  Figure  2.8  is  the 
corresponding  U  and  V  velocity  data  versus  the  time  from  release  for  partilce  #5.  The  variables 
figure  2.7  are  the  horizontal  and  vertical  distance  from  the  release  point  normalized  with  the  distance 
the  particle  could  fall.  A  small  bounce  is  indicated  in  the  figures  and  one  or  two  bounces  were 
typical  to  most  of  the  particles.  These  trajectory  and  velocity  history  data  are  similar  to  the  computed 
results  in  Part  III  of  this  report  but  the  computed  trajectories  consistently  impact  the  plate  upstream  of 
the  observed  impact  locations.  There  are  other  differences  between  the  observed  and  computed 
trajectories  which  are  commented  on  in  the  following  discussion. 

The  bounce  of  the  particles  after  impact  obviously  affects  the  disturbances  generated  at  the 
surface  and  possibly  turbulent  spot  formation.  However,  for  the  cases  near  the  critical  condition, 
there  is  no  definite  indication  that  spot  formation  is  affected.  For  the  polystyrene  particles,  a 
turbulent  spot  was  observed  with  particle  ^17  on  table  2.1  but  not  with  #18;  neither  particle  bounced 
detectably.  The  nylon  particle  #S  had  a  small  bounce  and  generated  a  spot.  Nylon  particle  #6  may 
have  bounced  but  this  was  small  enough  that  it  was  difficult  to  detect  and  no  spot  was  formed.  The 
two  delrin  particles  straddling  the  critical  condition,  #11  and  #12,  both  bounced  quite  noticeably. 
Thus,  all  combinations  of  bouncing  and  not  bouncing  subcritical  and  supercritical  particles  were 
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observed.  When  spots  did  form,  they  emerged  1  s  to  1 .33  s  after  impact  in  all  cases  and  there  is  no 
clear  indication  the  bounce  had  an  effect  on  spot  formation. 

Some  of  the  particle  velocity  histories  exhibited  oscillatory  behavior  that  was  not  predicted  in 
the  computations  in  Part  III.  Figures  2.9  and  2.10  are  velocity  histories  of  some  of  these  more 
unusual  cases,  particles  #7  and  ^  on  table  2.1.  The  pathlines  for  these  two  cases  are  the  delrin  data 
shown  on  figure  2.7  and  the  pathline  slope  is  increasingly  negative  as  impact  is  approached.  This 
trend  in  the  pathlines  indicates  that  the  U  component  of  velocity  is  decreasing  substantially  in  the  last 
stages  of  the  trajectory  relative  to  the  V  velocity  and  this  is  the  case  in  figures  2.9  and  2.10.  The 
velocity  history  for  sphere  ftl  shows  an  overshoot  of  the  freestream  velocity  of  about  20%  followed 
by  an  undershoot  of  U,  by  30%.  Sphere  #9  overshoots  U,  by  almost  40%  then  undershoots  it  by 
almost  50%  at  impact.  This  behavior  is  typical  to  the  delrin  particles.  There  were  similar  variations 
for  the  smallest  two  delrin  particles  which  overshoot  and  then  undershoot  U,  by  25%  and  30%  but 
then  relaxed  back  to  a  velocity  near  U«  at  impact. 

The  larger  diameter  nylon  sphere  data  were  similar  to  the  delrin  data  with  noticeable  velocity 
over  and  undershoots.  The  smaller  nylon  spheres  and  the  larger  polystyrene  spheres  were  similar  to 
particle  #5  data  shown  in  figures  2.7  and  2.8  with  the  U  velocity  data  consistently  exceeding  U.  by 
as  much  as  15%.  The  smaller  polystyrene  data  accelerated  to  a  velocity  within  10%  of  U.  with  both 
positive  and  negative  deviations  from  U.;  this  bdiavior  corresponds  with  the  estimated  accuracy  of 
the  measurement.  Thus,  there  is  a  consistent  trend  of  an  increasing  propensity  for  large  velocity 
oscillation  with  increasing  particle  diameter  and  increasing  specific  gravity.  The  smallest  and  lightest 
particles  produced  the  expected  results  which  includes  a  monotonic  trend  in  the  velocity  history  with 
no  deviations  of  the  streamwise  component  of  velocity  with  the  freestream  velocity  that  are  larger 
than  the  expected  error  as  the  particle  approaches  equilibrium  prior  to  impact. 
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The  vertical  velocities  all  rapidly  increase  to  near  the  impact  value  in  25%  or  less  of  the  time 
between  release  and  impact  for  all  diameters  and  materials.  Assuming  the  vertical  velocity  is  the 
impact  velocity  is  reasonably  accurate  for  the  purpose  of  discussion.  Also,  most  of  the  particles  have 
U  velocities  within  25%  of  U«  for  the  majority  of  their  trajectories.  Assuming  this  for  the  purpose 
of  a  general  discussion,  the  slip  velocity  Reynolds  numbers  typical  to  most  of  a  particle’s  trajectory 
can  be  estimated  from  impact  velocities.  The  smaller  polystyrene  particles  had  Reynolds  numbers 
based  on  the  two  component  slip  velocity  below  330.  Based  on  fixed  sphere  data  described  by 
Sakamoto  and  Haniu  (1990),  these  particles  are  shedding  hairpin  vortices  periodically  with  a  fixed 
orientation.  The  next  two  larger  polystyrene  particles,  #14  and  #15,  the  smaller  3  nylon  particles,  # 
4,  S,  and  6,  and  the  smallest  2  delrin  particles,  #1 1  and  12  had  slip  velocity  Reynolds  numbers  from 
370  to  840.  The  wake  may  not  be  turbulent  but  the  orientation  of  the  hairpin  vortices  shed  by  these 
spheres  is  likely  to  be  increasingly  irregular  with  increasing  Reynolds  number  in  these  cases.  The 
largest  polystyrene  particle,  the  larger  nylon  particles  and  most  of  the  delrin  particles  had  Reynolds 
numbers  near  1000  or  larger  and  a  turbulent  hairpin  vortex  shedding  process  with  an  irregular  and 
random  orientation  is  expected.  The  acceleration  of  the  spheres  in  the  present  experiment  involves 
both  the  speed  and  direction  and  this  may  change  the  character  of  the  wake  development  from  these 
fixed  particle  observations  of  Sakamoto  and  Haniu(1990)  but  it  is  expected  the  fixed  particle  behavior 
is  a  useful  guideline. 

The  unsteady  particle  wake  produces  unsteady  forces  on  the  particle  that  could  affect  its  path. 
It  is  expected  that  the  conditions  for  which  fixed  spheres  are  shedding  hairpin  vortices  regularly  with 
a  fixed  orientation  would  produce  the  least  scatter  in  the  pathline  and  velocity  history  data.  This 
apparently  occurs  for  the  smaller  and  lighter  particles  like  the  polystyrene  spheres.  Irregular  shedding 
with  random  orientations  and  finally,  turbulent  shedding  with  an  unstable  shear  layer  shed  from  the 
sphere  should  increase  the  spread  in  the  pathline  and  velocity  history  data. 
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The  computed  particle  trajectories  and  velocity  histories  in  Part  III  do  not  show  the  variations 
and  oscillations  observed  in  the  experiments.  The  model  used  in  Part  III  does  account  for  the 
transient  nature  of  the  particle  release  but  only  in  an  averaged  sense.  The  unsteady  character  of  the 
wake,  particularly  when  the  wake  is  expected  to  have  a  random  and  irregular  influence  on  the  drag 
force  on  the  particle,  is  not  incorporated  into  the  computations.  Consider  a  simple  example  that  the 
computed  model  would  fail  to  handle.  Spherical  particles  released  from  rest  in  a  water  filled 
isothermal  tank  with  no  flow  will  not  typically  strike  the  bottom  of  the  tank  at  a  point  exactly  below 
the  release  point  when  the  conditions  are  generally  similar  to  the  experiments;  transitional  and 
turbulent  wakes  with  specific  gravities  from  1.05  to  1.41.  The  mean  impact  position  on  the  tank 
bottom  for  a  large  number  of  trials  should  be  directly  below  the  particle  release  point  and  it  is  the 
mean  that  the  model  in  Part  III  predicts. 

The  potential  exists  for  velocity  overshoots  and  undershoots  if  the  unsteady  wake  can  have  a 
significant  affect  on  the  particle  path  and  there  is  an  appropriate  time  lag  for  the  wake  to  adjust  to  the 
changes  in  the  flow  conditions.  The  low  frequency  mode  Strouhal  number  for  the  hairpin  vortex 
shedding  frequency  of  a  fixed  sphere  in  a  uniform  flow  due  the  the  instability  of  the  averaged  wake  is 
about  0.2  for  700  <  Re^  <  3000  where  St  =  fdAJ,  see  Sakamoto  and  Haniu(1990).  Using  the  slip 
velocity  and  this  value  of  the  Strouhal  number,  shedding  frequencies  range  from  S  to  10  Hz  for  the 
nylon  and  delrin  particles  with  the  lower  values  for  the  larger  particles.  The  polystyrene  particle 
shedding  frequencies  are  less  than  3  Hz.  The  time  from  release  to  impact  spanned  a  range  from  0.6  s 
to  0.8  s  for  delrin  particles  to  2.0  s  to  2.3  s  for  the  polystyrene  particles.  The  various  particles  shed 
4  to  10  vortices  in  time  span  between  release  and  impact  based  on  these  estimates  from  steady  flow 
fixed  sphere  data.  Since  the  process  of  shedding  one  vortex  effects  the  development  of  the  next,  a 
wake  adjustment  time  to  the  changing  particle  velocity  of  two  shedding  cycles  may  be  reasonable.  In 
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the  typical  case,  the  particle  travels  30%  of  its  path  in  this  time.  This  adjustment  time  is  comparable 
to  the  time  span  of  the  swings  in  the  U  component  velocity  associated  with  the  over  and  undershoots. 

The  vertical  velocity  data  do  exhibit  irregularity  but  this  is  not  like  the  low  frequency 
oscillations  seen  in  the  streamwise  particle  velocity  histories.  The  vertical  velocity  changes  are 
generally  smaller  and  at  higher  frequency.  Abrupt  changes  in  the  acceleration  in  the  vertical  velocity, 
indicated  by  changes  in  the  slopes  of  the  velocity  histories,  see  figures  2.8,  2.9,  and  2.10,  are  typical. 
The  terminal  velocities  for  a  particle  in  a  still  fluid  listed  in  table  2.1  all  exceed  the  impact  velocity 
but  the  particles  appear  to  be  accelerating  very  slowly  in  the  last  half  of  the  trajectories  prior  to 
impact  in  all  cases.  This  almost  asymptotic  approach  to  the  terminal  velocity  is  similar  to  the  vertical 
velocity  histories  computed  in  Part  III. 

The  U  component  particle  velocity  is  influenced  directly  by  the  unsteady  fluid  dynamic  forces 
on  the  particle  but  the  vertical  velocity  is  influenced  directly  by  the  constant  gravitational  pull  and  as 
well  as  the  drag  forces.  The  relative  magnitudes  of  the  fluid  dynamic  and  gravitation  forces  change 
with  particle  diameter  and  material.  Consider  a  middle  size  particle,  d  =  S.Sl  mm,  fixed  in  a  flow 
with  the  freestream  velocity  prior  to  release,  Re^  =  847.  As  a  simple  approximation,  an  empirically 
based  drag  coefficient  can  be  used  to  estimate  the  drag  force  on  the  sphere  for  a  uniform  steady  flow. 
From  White(1974),  =  F/('ApU"A)  =  (24/Rii^  +  {6/(1  +  Re,’^)}  +  0.4  =  0.628  for  the  d  = 

5.51  mm  particle  fixed  in  the  freestream.  The  drag  force  is  F  =  1.53  x  KT  Newtons.  This  is  the 
largest  drag  force  expected  on  the  particle  and  estimates  conditions  at  release.  The  gravitational  force 
is  4.28  X  10■^  1.46  X  ICT,  and  3.51  X  l(f*  Newtons  for  polystyrene,  nylon,  and  delrin  panicles. 

The  gravitational  force  approximately  equals  or  exceeds  by  a  factor  of  2  the  largest  fluid  dynamic 
drag  force  expected  with  nylon  and  delrin  for  this  middle  size  case.  For  the  largest  panicles  the 
gravitational  forces  are  relatively  larger  and  for  the  smaller  panicles  it  is  relatively  smaller.  Given 
that  the  gravitational  force  is  able  to  accelerate  all  of  the  panicles  to  close  to  the  impact  velocity  in 
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only  25%  of  the  of  the  time  from  release  to  impact,  it  seems  reasonable  that  unsteady  fluid  dynamic 
forces,  which  can  be  comparable  in  magnitude  to  the  gravitional  force,  could  cause  the  over  and 
undershoots  of  the  freestream  velocity.  These  force  estimates  also  indicate  that  it  is  probably  the 
character  of  the  wake  and  not  the  ratio  of  the  gravitational  to  fluid  dynamic  drag  that  is  important  in 
the  present  cases.  The  particles  with  the  relatively  larger  gravitational  forces  have  the  most  irregular 
behavior  and  the  turbulent  wakes. 

Other  possible  explanations  for  the  U  velocity  history  oscillations  have  been  sought  simply 
because  the  observed  U  velocity  variations  were  not  expected.  Certainly  there  is  concern  that  errors 
in  the  procedure  could  be  a  contributing  factor.  In  light  of  the  fact  that  the  particles  with  low  Re^ 
behaved  as  expected  and  the  faster  higher  Re^  particles  behaved  unexpectedly,  this  seems  unlikely. 
The  displacement  between  video  frames  is  larger  for  the  faster  particles  and  the  relative  uncertainty  in 
estimation  of  this  displacement  decreases  as  a  result.  Repeated  trials  would  have  helped  to  define  the 
range  of  variation  which  is  not  known  based  on  the  current  experiments.  Thus,  it  is  not  possible  to 
judge  how  extreme  some  of  the  variations  observed  really  are.  Undetected  defects  in  the  particles, 
attached  small  bubbles  that  were  not  noticed,  or  the  influence  of  the  release  process  and  the  probe 
wake  are  other  possible  explanations  for  the  U  component  velocity  oscillations.  The  topic  of  pathline 
and  velocity  history  irregularity  for  falling  spheres  in  a  uniform  quiescent  freestream  requires  further 
study. 

2.5  Summary  and  Conclusions 

Flow  visualization  and  particle  tracking  have  been  used  to  study  turbulent  spot  generation  by 
spherical  particles  that  fall  freely  into  a  laminar  boundary  layer.  The  critical  Reynolds  number  based 
on  the  velocity  of  the  particle  at  impact  and  the  particle  diameter  is  approximately  700  to  800.  This 
is  in  agreement  with  results  on  a  laminar  flow  control  heated  body  by  Lauchle,  et  al(1986)  and 
supports,  for  values  of  I.85>  d/5*  >  0.5,  that  the  critical  Reynolds  number  is  constant.  Although 
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this  is  similar  to  the  critical  Reynolds  number  that  fixed  surface  mounted  spheres  generate  turbulent 
wedges,  the  processes  involved  do  not  appear  to  be  similar.  In  all  cases,  the  turbulent  spots  emerge 
in  the  inunediate  vincinity  of  the  convecting  particles  and  it  appears  that  this  is  a  result  of 
disturbances  created  by  the  particle  near  or  below  y/5  =  0.2.  The  larger  particles  convected 
downstream  on  the  surface  with  the  turbulent  spots  they  generated  but  the  smaller  particles  rolled  on 
the  surface  of  the  plate  at  lower  velocity  than  the  spots  they  generated.  In  no  case  were  particles 
observed  to  generate  more  than  one  turbulent  spot  or  to  create  a  turbulent  wedge  after  impact.  The 
spots  begin  to  emerge  within  a  time  from  impact  of  4S  <  tU./5*  <  60,  approximately,  and  the  spots 
emerge  with  approximately  25  particle  diameters  the  impact  location.  Turbulent  wedges  form  in  the 
wake  of  a  fixed  sphere  SO  to  75  particle  diameters  downstream  of  the  particles  at  the  critical 
condition.  In  some  of  the  convecting  particle  cases  with  less  complicated  and  smaller  appearing 
disturbances  at  impact,  spanwise  oscillations  are  observed  after  impact  and  prior  to  the  emergence  of 
a  turbulent  spot.  These  oscillations  may  indicate  a  process  of  selective  amplification  but  the  short 
time  span  between  the  particle  impact  and  the  emergence  of  a  turbulent  spot  suggests  this 
amplificattion  stage  is  nonlinear  and  a  bypass  transition  process  is  active. 
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(b) 

Fig.  2.1  Flow  visualization  of  d  =  5.56  mm  diameter  polystyrene  particle  generated  events, 
(a)  t  =  0  s,  (b)  t  =  0.66  s 
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(d) 

Fig.  2.1  Flow  visualization  of  d  =  5.56  mm  diameter  polystyrene  particle  generated  events. 
(c)t  =  1.33  s,  (d)t  =  2.0  s 
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Fig.  2.2  Impact  of  a  subcritical  d  =  3.97  mm  diameter  polystyrene  sphere,  (a)  t  =0 
s,  (b)  t  =  Is,  (c)  t  =  2s. 


60 


Fig.  2.3  Events  generated  by  impact  of  d  =  6.35  mm  diameter  polystyrene  particle 
at  (a)  t  =  0  s,  (b)  t  =  0.66  s,  and  (c)  t  =  1.33  s. 
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Fig.  2.S  Two  component  impact  velocity  magnitude,  U^,  versus  Reynolds  number. 
Solid  symbols  are  conditions  that  did  not  generate  a  turbulent  spot. 
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Fig.  2.6  Estimated  critical  conditions  for  turbulent  spot  generation  on  a  heated 
laminar  flow  control  body  at  three  heating  conditions  ft^om  Lauchle,  et  al.(1986). 
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Fig.  2.7  Trajectories  for  particles  #S,  7,  and  9  on  table  2.1.  h  is  the  distance  the 
particle  can  fall  from  the  release  point. 
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Fig.  2.8  Velocity  history  for  particle  #5  on  table  2.1,  nylon  with  d  =  3.94  mm. 
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1.9  Velocity  history  of  particle  #7  on  table  2.1,  delrin  d  =  7.89  mm. 
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Table  2.1  Impact  conditions  for  all  cases  considered. 


ijjjlll 

d  (mtn) 

Material 

d/Sit 

Impact 

uju. 

V./U. 

Re^  at 
impact 

RCta,  at 
impact 

Estimated 

terminal 

velocity, 

V./U. 

Spot  on 
impact 

1 

7.89 

Nylon 

1.03 

0.56 

-0.96 

1170 

1350 

-1.74 

yes 

2 

6.31 

Nylon 

0.80 

0.93 

-1.19 

1160 

1470 

-1.58 

yes 

3 

5.51 

Nylon 

0.69 

1.10 

-0.99 

837 

1250 

-1.55 

yes 

4 

4.72 

Nylon 

0.58 

0.94 

-1.01 

736 

1060 

-1.36 

yes 

5 

3.94 

Nylon 

0.48 

0.88 

-0.82 

496 

726 

-1.29 

yes 

6 

3.12 

Nylon 

0.37 

0.82 

-0.77 

371 

542 

-1.14 

no 

7 

7.89 

Delrin 

1.16 

0.73 

-2.10 

2550 

2700 

-2.97 

yes 

8 

6.31 

0.86 

1.12 

-1.96 

1900 

2190 

-2.58 

yes 

9 

5.51 

0.77 

0.49 

-1.67 

1410 

1470 

-2.38 

yes 

10 

4.72 

0.64 

0.92 

-1.61 

1170 

1350 

-2.29 

yes 

11 

3.94 

0.54 

0.93 

-1.36 

824 

998 

-1.95 

yes 

12 

3.12 

Delrin 

0.42 

0.95 

-1.14 

547 

711 

-1.59 

no 

13 

7.94 

Polystyrene 

0.92 

1.08 

-0.81 

987 

1640 

-0.98 

yes 

14 

7.14 

Polystyrene 

0.82 

0.98 

-0.70 

769 

1320 

-0.92 

yes 

15 

6.35 

Polystyrene 

0.71 

0.99 

-0.66 

646 

1160 

-0.85 

yes 

16 

5.56 

Polystyrene 

0.62 

1.03 

-0.63 

537 

1030 

-0.78 

yes 

17 

4.76 

Polystyrene 

0.51 

0.99 

-0.43 

317 

791 

-0.71 

yes 

18 

3.97 

Polystyrene 

0.42 

1.04 

-0.54 

330 

787 

-0.63 

no 

19 

3.18 

Polystyrene 

— 

0.96 

-0.48 

235 

523 

-0.54 

no 
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PARTin 

Particle  Trajectory  and  Transition  Analysis 
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1.  Introduction 


The  analytical  and  computational  approach  taken  in  this  research  progiain  lias 
been  guided  b}'  the  results  of  the  parallel  experiments.  In  the  originall}-  pioposed 
technical  approach  it  was  planned  to  determine  the  trajectory  of  small  freely- 
suspended  particles.  The  wave-packets  generated  by  these  particles  were  to  be 
simulated  by  moving  point  disturbances  using  the  method  developed  by  Brener 
and  Haritonidis  (1990).  In  this  method,  the  linearized  equations  of  motioned 
are  transformed  in  space  and  marched  in  time.  The  physical  flow  field  at  any 
time  is  recovered  by  an  inverse  transformation.  This  approach  avoids  some  of 
the  difficulties  of  a  normal  mode  or  eigenfunction  expansion  technique  in  the 
description  of  the  initial  conditions.  In  the  proposed  calculations,  both  foi  the 
particle  trajectory  and  wave-packet  development,  the  basic  flow  was  taken  to  Ire 
the  Blasius  boundary  layer;  however,  the  accompanying  experiments  indicated 
that  the  mean  velocity  defect,  found  in  the  particle's  wake,  was  significant  and 
played  an  important  role  in  the  determination  of  the  initial  disturbance  and  its 
subsequent  development.  Thus,  it  was  decided  to  include  this  effect  in  the  analysis 
and  computations. 

To  this  end,  the  stability  of  the  wake  of  a  body  suspended  in  a  boundary 
layer  has  been  examined.  This  has  involved  several  steps.  First,  the  stability 
characteristics  of  a  two-dimensional  wake  in  a  Blasius  boundary  layer  have  been 
examined.  The  role  of  the  wake  position,  width,  and  velocity  defect  and  the  rela¬ 
tive  importance  of  viscous  and  inviscid  or  dynamic  instability  have  been  studied. 
As  a  second  step,  the  same  problem  is  being  solved  with  a  perturbation  analysis. 
In  this  analysis,  the  small  parameter  is  the  strength  of  the  mean  shear  that  sur¬ 
rounds  the  wake.  This  analysis  leads  to  the  third  problem:  the  stability  of  the 
wake  of  a  sphere  embedded  in  a  shear  flow.  Through  the  use  of  a  perturbation 
analysis  this  problem  may  be  solved  to  leading  order  as  an  axi.symmelric  wake 
in  a  uniform  flow.  The  effects  of  the  shear  flow'  are  included  at  higher  order  as 
corrections  to  the  axisymmetric  problem.  This  approach  avoids  the  need  to  solve 
a  non-separable  boundary  problem  for  the  stability  of  an  axisymmetric  wake  in  a 
plane  shear  flow.  The  results  of  the  final  part  of  this  analysis  provide  the  initial 
conditions  for  the  development  of  a  wave-packet  generated  by  a  ])ariicle  suspended 
in  or  traversing  a  boundary  layer.  In  addition,  the  problem  is  related  closely  to 
the  experimental  configuration  and  direct  comparisons  betw'een  the  predictions 
and  experiments  will  be  possible. 

Thus  the  analysis  and  computations  include  the  following  component  prob- 
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lems; 


•  the  calculation  of  the  trajectory  of  freely-suspended  particles  in  a  laininai 
boundary  layer. 

•  the  prediction  of  the  development  of  wave  packets  in  a  laminar  boundai  \ 
layer  generated  by  disturbances  at  different  locations. 

•  the  calculation  of  the  stability  characteristics  of  a  two-dimensional  wake  em¬ 
bedded  in  a  laminar  boundary  layer  with  a  full  numerical  and  a  ix-rt  iii  l>aiion 
analysis. 

•  The  prediction  of  the  stability  of  the  wake  of  a  sphere  embedded  in  laminar 
boundary  layer,  and 

•  the  development  of  wave  packets  generated  by  finite-size  particles  susjrended 
in  a  laminar  boundary  layer. 

The  first  two  items  above  have  been  completed  as  has  the  numerical  solution 
of  the  third  item.  These  are  described  in  this  report.  The  remaining  items  aie 
being  completed  as  parts  of  the  doctoral  dissertation  of  ms.  .Anupa  Bajwa.  The 
completion  of  this  work  is  expected  before  December  1993.  A  copy  of  the  thesis 
and  any  companion  publications  will  be  provided  to  the  sponsor  at  that  lime. 
Thus  the  present  report,  though  a  final  report  for  formal  reporting  purposes,  does 
not  represent  the  completion  of  our  activities  in  this  subject. 

In  the  next  section  the  particle  trajectory  calculations  are  presented.  This 
is  followed  by  a  description  of  the  wave-packet  development  from  infinitesimally 
small  particles  suspended  in  a  laminar  boundary  layer.  The  final  section  describes 
the  stability  characteristics  of  a  wake  embedded  in  a  laminar  boundary  layer. 

2.  Particle  Trajectory  Calculations 

.A  particle  released  into  a  fluid  is  acted  on  by  various  bod\’  and  surface  forces 
which  result  in  the  motion  of  the  particle  through  the  fluid.  The  trajector_\-  of  such 
particles  is  to  be  computed  to  determine  the  locations  and  velocities  of  spherical 
particles  of  various  diameters,  specific  gravities  and  release  velocities. 
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2.1.  Second-order  boundary  layer  theory 

A  freely  suspended  particle  in  a  fluid  is  known  to  affect  boundary  layer  transition. 
To  determine  the  trajectory  of  such  a  particle,  the  basic  flow-field  (n.f  )  must  l)e 
computed  first.  For  a  uniform  freestream  the  first-order  boundary  layer  llieor\ 
gives  a  zero  vertical  velocity  outside  the  boundary  layer.  Hence  second-order 
boundary  layer  effects  have  to  be  taken  into  account,  using  the  theory  of  matched 
asymptotic  expansions  (Van  Dyke.  1975). 

For  a  flat  plate,  the  flow  due  to  displacement  thickness  is  given  by  the  com¬ 
posite  stream  function  v(x.y:R)  : 


I  -  f  ( '^y\  _  -^1 

Vcon,pos.te-  ^ ^  ^ 


+  i  !/  '. 


where  R  is  the  Reynolds  number  of  the  flow,  and  3?  denotes  the  real  part 
complex  function.  Therefore  the  velocity  components  are  given  bv 
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Hence  the  local  shear  components  can  be  analytically  obtained  from. 
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From  these  equations,  the  convective  derivatives  for  U  and  V  are  computed  as 


DU  ^  ^ 

Dt  dx  dy 


DV  ^  ^ 

Dt  dx  dx 


These  quantities  are  required  to  compute  the  particle  trajectory. 
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2.2.  Trajectory  calculations 

Once  the  basic  velocity  field  is  computed,  the  particle  trajectory  can  be  deteniiiut'd 
as  follows: 

For  a  spherical  particle  p  of  radius  a*,  the  motion  through  a  viscous  fluid  f  is 
governed  by  the  equation  (Hoffman,  198S) 


df  '  r 

.dV'  dv: 


where  =  particle  volume  =  .4*  =  particle  surface  area  = 

The  various  terms  in  this  equation  represent: 


LHS  :  Mass  of  the  particle  multiplied  by  its  acceleration 
RHSj:  Buoyant  force  on  the  particle 

RHSj;  Force  due  to  the  pressure  gradient  in  the  fluid  surrounding  the  particle 

RHS3:  Drag  force  on  the  particle 

RHS4:  Lift  force  on  the  particle 

RHS5:  Force  required  to  accelerate  the  added  mass 

RHSe:  Basset  force  on  the  particle,  which  dominates  during  tlie  initial  accelera¬ 
tion 


The  variables  may  be  non-dimensionalized  with  respect  to  an  arbitrary  length 
scale  L'  and  the  freestream  velocity  .  Then,  the  particle  trajectory  equation 
becomes. 


2^  Pf  dt 


(9) 


where  Fr  is  the  Froude  number  =  jg’L',  R  is  the  Reynolds  number  = 
U^L’ /u’  ,  and  Is  is  the  Basset  integral  =  /o  l\/f  —  ~  ■ 
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The  system  of  equations  for  the  trajectory  is  completed  b\  specifyin  g  tin- 
position  vector  fp,  where. 


_  i7 

dt  ~  " 


(10) 


For  flow  along  a  horizontal  flat  plate  the  component  equations  are. 
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For  a  given  Reynolds  number,  radius,  specific  gravity,  initial  location  and  veloc¬ 
ity  of  the  particle,  these  equations  may  be  integrated  in  time  to  determine  its 
trajectory. 


2.3.  Lift  and  Drag  Terms 

It  can  be  seen  in  Figure  1  that,  for  a  typical  case  of  a  freely  released  particle,  the 
lift  and  drag  forces  have  parallel  y-components  and  anti-parallel  j-components. 
The  drag  coefBcient  depends  on  the  particle  Reynolds  number  Ftp  =  2a:lqR  as 
follows: 

logioCc  =  1.6435-1.1242(log,oRp)+0.158S(log,oRp)^  where  260  <  Rp  <  1500 

(IS) 
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Figure  1:  Lift  and  drag  force  vectors  on  a  particle 


This  is  the  recommended  drag  correlation  based  on  the  standard  drag  curve  (Clifi. 
Grace  and  Weber,  1978).  In  the  range  of  Rp  expected  in  the  present  exijcrimenis. 
the  lift  coefficient  depends  only  on  the  local  fluid  shear  (Dandy  and  Dwyer.  1990). 
Since  dV/dx  is  very  small  throughout  the  flow  regime,  we  set  CT,  =  0  and 

=  0.30o*{dU/dy)p  for  Rp>50  where  {dU/dy)^^  =  dUjdy*a/^(j  (19) 


2.4.  The  Basset  Term 


The  final  term  on  the  right  hand  side  of  the  equations  of  motion  is  the  Basset 
history  integral  in  which  the  past  acceleration  of  the  particle  is  included.  Tlic 
integrand  of  the  Basset  term  is  singular  at  the  upper  limit  t  =  t.  However,  liie 
Basset  force  can  be  included  in  the  calculations  using  a  method  which  eliminates 
the  need  to  evaluate  the  integral  explicitly.  A  change  of  variables,  as  suggested 
by  Hoffman  (1988),  allows  the  trajectory  equation  to  be  written  in  a  simple  form 
as  shown: 

As  =  Bn  +  CnIbn  time  is  (20) 


where 


y  di  dt  j 


(21 
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Thus  1b  can  be  eliminated  and  the  particle  acceleration  at  is  can  be  found  as 
follows: 

/  —  \  y  _ 

DV 


Dt 
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where 


As  = 


Bs  ~  C y/tsh 
1  +  C\/tyv  A(p2 


(27) 


We  have  used  this  trajectory  formulation  to  conduct  a  parametric  study  of 
particle  trajectories  for  a  range  of  particle  properties  and  release  points  to  be 
considered  in  the  experiments. 

The  (non-dimensional)  parameters  considered  are 


Reynolds  number  R  =  250,000 

Particle  radius  a  =  0.001,0.002,0.003,0.004 

Density  ratio  Pp/ pj  =  1.05, 1.15. 1.25, 1.35 

Release  location  Xp(l)  =0.1.  yp(l)  =  0.02,0.1 

Release  velocity  Up{l)  =  0.0, 0.5, 0.9, 1.2,  lp(l)  =  0.0 
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In  the  calculation  of  the  trajectories  for  the  various  cases,  it  is  oljserwd  that 
the  buoyancy  force  atid  drag  dominate  the  motion  throughout  the  flow  fudd.  and 
the  lift  force  is  significant  only  in  the  boundary  layer,  where  it  reaches  a  niaxiinuni 
value  of  S9c  of  the  drag.  .A,s  expected  from  the  analysis  of  Saffman  ( 196')).  lot  a 
p.art’icle  moving  slower  than  the  fluid  over  a  stationary  wall  tlie  lift  lends  to  pu^ii 
the  particle  away  from  the  wall.  In  the  case  of  plate  mo\  ing  through  a  siaiiunai  \ 
fluid,  the  situation  could  be  reversed  depending  on  the  relative  motion  between 
the  particle  and  the  fluid  (Fig.  2).  This  is  best  understood  when  viewed  from 
the  frame  of  reference  of  the  particle.  The  added  mass  term  and  the  Basset  force 
are  found  to  not  be  significant  in  the  calculations,  since  the  density  ratio  does  not 
tend  to  zero  for  the  case  of  solid  spheres  moving  through  a  fluid.  I  hesc  fijicv- 
would  be  quite  sigiuficant.  for  example,  in  the  calculation  of  a  bubble  trajeclm  w 
Hence,  for  the  particle  trajectories  (Figs.3-(>).  the  Basset  term  is  omitted  from 
the  calculations  in  order  to  simplify  them. 

Fig.  3  shows  the  variaiton  of  the  trajectories  of  freely  released  partic  les  for 
different  particle  densities. Fig.  4  shows  the  varaition  with  particle  radius.  Figures 
5  and  6  show  the  corresponding  time  variation  of  the  particles'  vertical  velocity. 
As  expected,  the  smaller  and  lighter  particles  travel  through  a  greater  horizontal 
distance  and  impact  the  plate  further  downstream  than  the  bigge-r  and  heavic! 
ones.  Fig.  5.  shows  that  the  lighter  particle  reaches  a  smaller  terminal  velocity  and 
that  it  ta.kes  longer  to  reach  it.  implying  that  its  acceleration  vertically  downwards 
is  of  a  smaller  magnitude.  This  is  because  Ts  bod}’  force  is  smaller.  When  the 
particle  radius  is  varied,  as  in  Fig.G.  it  is  observ'ed  that  the  smaller  panicle  reaches 
a  lower  terminal  velocity  and  that  it  takes  less  time  to  do  so  than  the  biggei 
particles.  This  means  that  its  drag  force  increases  rapidly  in  time  to  iraiance  its 
weight.  This  is  because  the  drag  coefficient  Co  >5  higher  for  smaller  particles, 
initially,  and  it  increases  in  time  as  Rj,  decreases. 

3.  Development  of  wave  packets  in  a  laminar  boundary 
layer 

.A  particle  disturbs  the  flow  as  it  moves  along  its  trajectory.  So  the  next  stej)  is  to 
compute  the  wave  packet  development  from  this  particle  point-source.  To  do  this 
a  linear,  inviscid  stability  analysis  is  used,  for  a  slationai}’  particle  to  start  n  iih 
There  are  two  different  ways  to  approach  this  problerri  —  to  follow  theex  ohition  ol 
the  disturbance  in  time  or  to  compute  the  normal  modes  in  the  fre(|nenc\  domain 
and  superimpose  them. 
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a)  Flow  over  a  stationary  wall 


b)Motion  of  wall  through  a  still  fluid 


Laboratory  frame  of  reference 


Particle  frame  of  reference 


Figure  2:  Lift  force  on  a  particle  due  to  shear 
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Figure  4;  Trajectories  of  particles  with  different  sizes  released  from  rest 
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Figure  5:  Time  variation  of  the  vertical  velocity  of  particles  with  differenl  specific 
gravities 
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A  wave  packet  generated  by  an  ini|)ulse  disturbance  causes  usc  illaiion^  in  lli< 
fiow  as  it  spreads  out  while  travelling  downstream  ((lastei  ami  Ciraiil.  A)7'i; 
This  initial  disturbance  generates  all  possible  modes  and  the  wa\(  packet  lutni- 
through  selective  amplification  and  interference  of  the  most  unstable  waxes.  Thi^ 
phenomenon  can  be  simulated  numerically  using  a  linear  stability  analysis  (Ciastei . 
1975).  The  wave  packet  is  reproduced  computationally  by  linearly  combining 
spatially  unstable  modes  and  summing  over  all  wavenumbers  and  frec|uencies. 

To  study  the  initial  value  problem  for  any  general  disturbance  in  a  viscous 
flow,  the  Orr-Sommerfeld  equation  must  be  solved.  The  eigcnsolutions  to  this 
equation  are  discrete  modes  which  form  a  complete  system  for  bounded  fiow.s. 
However,  for  boundary  layer  flows  with  a  semi-infinite  domain.  theK-  is  also  a 
viscous  continuous  spectrum  (Gustavsson.  197S).  .As  the  Reynolds  number  in¬ 
creases.  one  mode  approaches  the  inviscid  Rayleigh  mode  while  the  higher  modes 
tend  towards  a  neutrallx’  stable  inviscid  continuum.  For  three-dimensional  dis¬ 
turbances  in  a  piecewise  linear  boundary  layer  the  vertical  velocitx'  has  two  parts 
-  one  is  a  dispersive  part  which  comes  from  the  Rayleigh  equation  solution,  and 
the  other  is  a  convected  part  from  inviscid  continuous  spectrum.  .Also,  there  is  a 
permanent  scar  caused  by  the  fact  that  the  streamwise  velocity  does  not  vanish 
at  large  times.  It  is  advected  at  the  local  mean  velocity.  .Most  three-dimensional 
disturbances  show  a  linear  increase  in  energy  with  time.  This  bcdiavioui  i--  dif¬ 
ferent  from  the  one  predicted  by  two-dimensional  stability  analysis  that  transicni 
modes  decay. 

Time-evolution  of  a  small  disturbance  can  be  computed  using  Fourier  trans¬ 
form  methods.  Results  show  that  transient  modes  do  not  decay  as  rapidly  as 
dispersive  modes  (Henningson,  1988).  This  method  is  extended  to  the  Blasius 
boundary  layer  by  Breuer  and  Haritonidis  (1990).  Both  numerical  and  experi¬ 
mental  results  indicate  the  formation  of  linearly  unstable  wave  packets,  and  of 
dominant  transients  which,  for  weak  initial  disturbances,  eventually  decay  due  u> 
viscous  effects.  For  weak  disturbances  a  linear  stability  analysis  is  applicable  and 
this  approach  is  discussed  in  detail  in  the  next  section.  For  strong  distm  !)anr<'s 
direct  numerical  simulation  is  necessary.  For  these  high-amplitude  disturbance^ 
the  transients  do  not  decay.  Non-linear  effects  lead  to  the  direct  Ijreakdown  of  tlie 
disturbance  to  a  turbulent  spot,  bypassing  the  wave  packet  stage. 

For  low-amplitude  initial  disturbances  in  a  flat-plate  boundary  la\  er  at  moder¬ 
ate  Reynolds  numbers,  the  route  to  transition  is  through  a  wave  packet  breakdown 
which  occurs  over  a  broad  spatial  domain  (Tso.  Chang  and  Blackwelder.  1990). 
Since  the  grox.ing  wave  packet  moves  significantly  slower  than  the  diminishing 
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transient,  hump-like  disturbance,  there  is  little  interaction  between  them  and  oni\ 
the  wave  packets  are  responsible  for  transition.  W  ith  strong  initial  disiui  Ijancc^ 
however,  the  transition  mechanism  is  via  the  growth  of  the  Ir'insuitlt^  into  a  lui- 
bulent  spot.  We  will  consider  only  w'eak  disturbances  in  modeling  the  effect  of 
the  particles. 

4.  Numerical  simulation  of  wave  packets 

Caster's  theoretical  linear  model  successfully  simulates  his  experimentally-oljserved 
wave  packet  excited  by  a  localized  pulse  at  the  wall.  While  Caster‘«  simple  lin¬ 
ear  model  works  well  for  calculations  at  a  fixed  height  above  the  plate,  it  cannot 
account  for  a  movement  of  the  source  in  the  vertical  direction.  A  recent  method 
developed  by  Breuer  and  Haritonidis  (1990)  employs  a  time-marching  teclmi<(ne 
to  compute  the  development  of  wave  packets  from  small-amplitude  disLurl)ances. 
This  approach  does  capture  the  dependence  of  the  wave  packet  on  the  y-locaiioii 
of  the  source  of  the  disturbance.  This  is  exactly  what  is  required  to  study  ilie 
effect  of  a  particle  disturbing  a  boundary  layer  as  it  moves  through  it.  We  use 
Breuer’s  analysis  in  the  development  of  our  computational  scheme.  To  verify  the 
code  we  compare  results  for  an  initial  disturbance  consisting  of  a  pair  of  countei  - 
rotating  vortices.  This  initial  disturbance  can  then  be  modified  to  model  the 
particle  moving  across  the  boundary  layer. 

The  equation  governing  linear,  inviscid,  small  disturbances  lor  the  vertical 
velocity  is  the  Rayleigh  equation 

(■^  -j-  iaU){^  ~  P)v  -  iqU"  V  =  0  (2S) 

dt  dy^ 

where  4-  0'^.  Here  a,0  are  the  streamwise  and  spanwise  wavenumbers 

respectively,  v  is  the  transformed  vertical  velocity  perturbation  obtained  by  taking 
a  Fourier  transform  in  x  and  z  of  the  physical  fluctuation  quantity  v: 

v{y,t;a,0)  =  f  [  vix,y,2,t)expl-iiQX  +  3z)]d.r  dz  (29) 

The  laminar  boundary  layer  profile  is  used  for  the  basic  flow’  U(y).  The  velocity 
fluctuations  perpendicular  to  and  along  the  wavefront,  and  |  re.speclix’ely.  arc 
related  to  v  bv 

(—  +  iaU)ui  =  +  laO’fi'y  i'-iO) 
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—  +  'of  )«’i  =  -Tt  '■ 
dt  k 

The  two  physical  velocity  perturbations  are  retrived  from  these  by 

i  =  (qui  +  dii'i )/ A- 
w  =  (Siij  —  ou'i  )/A' 

The  Rayleigh  equation  is  rewritten  as  a  system  of  two  equations 

d  f 


dt 


=  m(f  ."T  -  f  V 


9  ,  9 

V^k  =  i---k^)r 
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These  equations  are  marched  forward  in  time  using  a  Crank-Xicholsuu  scheme 
and  a  second-order  accurate  finite  difference  for  the  y-derivative.  for  which  the 
discretized  equations  are 

.2  ,-,n+l  _  „2 -n  ^  ^  -n)  _  {yV  (y  2 -n+1  ^  ^  (3(.  ^ 


V  -  V  Vj  = 


2.-,n  _  ^'j+1  _  jL2-n 


v‘t-; 


Ay2 


where  n  refers  to  the  current  time-step  and  j  refers  to  the  \ei  lical  grid.  In  \ecior 
form,  at  time  step  n  a  matrix  equation  is  obtained. 

In  vector  form,  at  time  step  n  a  matrix  equation  is  obtained 


where 


(I  -f  fR)Dv"+'  =  (I  -  tR)Dv"  +  zS(v"+‘  +  v"; 

v"  =  lf,T  ;  =  1.2. 

R=^I 


J 


s  = 


2 

oAtfT. 


(37) 

(3S) 

(39) 

(301 


I  is  the  identity  matrix  and  D  is  the  tridiagonal  matrix  associated  with  the  Lapla- 
cian  derivative  in  the  normal  direction. 
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Therefore  the  explicit  form  for  the  velocity  at  the  (t!  +  1  )ih  liine-steij 


=  (D  +  2(RD  -  S))-*(D  -  ?(RD  -  S))v" 


Similarly  for  i/'i  the  matrix  equation  is 

(I  +  2R)w"+i  =  (I  -  iR)w"  -  T( +  v’‘ ) 


(42) 


where 


T  = 


^i3U' 

2h 


(42) 


This  is  also  integrated  forward  in  time  using  the  current  values  for  v"  and  . 
The  quantity  ui  is  obtained  simply  from 


u,  = 


As  a  test  case  for  the  our  code,  we  use  the  initial  conditions  of  Brener  et 
al  which  simulate  a  localized,  impulsive  disturbance  in  the  form  of  two  pairs 
of  counter-rotating  longitudinal  vortices,  represented  b\’  the  (two-dimensional) 
stream  function 

'I' =  zf^^exp(— (45) 


where  x,y,z  are  the  three  Cartesian  coordinates  scaled  by  appropriate  scaling 
factors 

x  =  x/li\  y  =  y/ly-,  S=zlL  (-1C) 

From  this  stream  function,  the  initial  disturbance  velocities  are  computed  as 


V  =  w  =  (4  I ) 

dz  dy 

At  f  =  0  the  center  of  the  disturbance  is  located  at  i  =  0,r  =  0.  Since  all  the 
matrices  are  tridiagonal,  the  system  of  two  matrix  equations  is  easily  integrated 
forward  in  time.  Furthermore,  for  fixed  a,  0  and  At  the  matrices  need  to  be 
evaluted  and  combined  only  once.  Hence  time-marching  for  v  at  the  next  time- 
step  reduces  to  a  simple  multiplication  of  a  matrix  by  the  vector  v  at  the  current 
time  step.  This  speeds  up  the  computations  for  each  pair  of  [o.B).  64  modes  are 
used  for  a  and  32  for  0.  These  are  sufficient  to  produce  smooth  contour  plots 
for  the  various  velocity  fluctuations,  obtained  at  any  desired  time  i  l^y  taking 
the  (inverse)  fast  Fourier  transform  of  the  transformed  velocities  to  convert  the 
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Figure  7:  Contours  of  initial  vertical  velocity  in  the  ly-plane  at  c  =  0.  Conlour 
spacing  is  .1  of  the  maximumu  initial  t>-perturbation 

fluctuations  back  to  physical  space.  This  is  done  at  each  y-localioii  lo  get  coin|>l(*U' 
picture  across  the  boundary  layer.  Considerable  effort  has  been  made  to  use 
optimal  two-dimensional  fast  Fourier  transform  routines  in  the  comple.x  domain 
and  efficient  matrix  inversion  routines,  since  these  are  used  at  each  y-locat  ion  in 
the  boundary  layer.  The  routines  are  adapted  from  Press  et  al  (19SG). 

The  evolution  of  the  vertical  velocity  disturbance  is  shown  in  Figs.  7-11.  The 
initial  disturbance  of  two  counter-rotating  streamwise  vortices  develops  into  one 
having  a  wave  packet  part  and  a  transient  part.  The  wave  packet  nature  is  seen 
in  the  X2  plane  and  the  transient  part  is  clearly  seen  in  the  xy  plane  where  the 
patch  of  low-speed  fluid  pulls  ahead  of  the  rest  of  the  disturbance. 

Thus  wave  packet  development  can  be  computed  using  linear  stabiiit\-  theorx. 
A  simple,  computationally-efficient  time-marching  scheme  is  sufficient  to  study 
evolution  of  fluctuation  velocities  in  the  spectral  domain.  Computations  for  weak 
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Figure  8:  Contours  of  the  vertical  velocity  in  the  xt/-plane  at  r  =  0  at  time  t 
Contour  spacing  is  .1  of  the  maximumu  initial  u-perturbation 
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Figure  9:  Contours  of  the  vertical  velocil 
t  =  25.  Contour  spacing  is  .1  of  the  maxi 
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Figure  10:  Contours  of  the  vertical  velocity  in  the  xy-plane  at  r  =  0  at  time 
t  =  50.  Contour  spacing  is  .1  of  the  maximumu  initial  i>-perturbalion 


initial  disturbances  can  be  extended  to  include  moving  source^,  mh  li  a^  pailitlc' 
in  a  boundary  layer.  For  a  moving  particle,  the  wave  packets  generated  at  each 
j/-location  can  be  computed  individually  and  then  sujrerposed  linearly  to  gel  the 
entire  flow  field,  which  is  possibly  a  turbulent  spot.  Work  is  in  progress  on  liuw 
exactly  to  model  the  particle  as  an  initial  distur'oance.  The  m<»del  is  based  mu^llx 
on  experimental  results  available  in  the  literature. 

4.1.  Initial  disturbance  due  to  a  particle 

To  model  a  freely-suspended  particle  as  an  initial  disturbance,  we  must  idenlib 
the  flow  pattern  immediately  downstream  of  the  particle  in  the  Ijoundarv  lavci 
flow.  We  model  a  spherical  particle  to  start  with.  \'ortices  slied  from  s|)licres  arc 
basically  of  two  types:  \ortex  tubes  and  large-scale  hairi)in  \oitici's.  .\  lc(l•nl 
experimental  study  by  Sakamoto  and  Haniu  (1990)  concludes  ihal  a  hich  type  ol 
vortex  is  shed  depends  strongly  on  the  Reynolds  number  of  the  (low.  For  otii 
range  of  Reynolds  numbers  (based  on  particle  diameter  and  freest  ream  \  eloc  ii>') 
the  primary  vortex  pattern  is  that  of  hairpin-like  vortices  of  known  Sirouhal  num¬ 
ber.  These  are  caused  by  the  progressive  wave  motion  of  llie  wake  and  ate  shed 
periodically.  The  amplitude  is  irregular  and  the  direction  of  shedding  oscillates  in¬ 
termittently.  This  may  be  due  to  an  imbalance  in  the  suppl  v.  slot  age  a/icl  cmis.Mm/ 
of  energy  within  the  vortex  formation  region.  The  length  scale  of  llu'sc  vortices 
can  be  determined  based  on  flow  velocity  and  shedding  frequency.  For  Reynolds 
numbers  higher  than  700  the  hairpin  vortices  become  turbulent  and.  in  addition 
to  these,  small  vortex  loops  are  shed.  Some  of  the  loops  move  downstream  and 
intertwine  with  the  hairpin  vortices. 

To  start  with,  we  will  model  the  hairpin  vortices  behind  thesphei  e  as  a  couiikm  - 
rotating  streamwise  vortex  pair.  For  a  particle  tethered  to  the  wall,  these  vort ice's 
are  generated  periodically  at  the  same  location  in  space.  So,  in  our  linear  stability 
calculation  of  the  wave  packet  we  can  simply  superpose  the  fluctuations  due  to 
several  initial  disturbances  staggered  in  time.  The  calculalioms  need  to  be  pei- 
formed  only  once,  and  thereafter  we  can  simply  shift  the  time-origin  to  include  t  he 
periodically-shed  vortices.  The  strength  of  the  initial  disturbance  changes  in  time. 
However,  the  vertical  velocity  fluctuation  can  be  scaled  out  in  the  spectral  domain. 
So  it  is  necessary  to  compute  the  evolution  of  only  a  known  strength  and  a  unit 
strength.  Then  for  any  other  values  of  initial  disturbance,  the  iime-e\olin ion  of 
the  fluctuation  can  be  obtained  by  a  simple  scaling. 

For  a  moving  particle,  the  location  of  the  initial  disturbance  changes  in  time. 
Hence  the  calculations  have  to  be  performed  using  each  vortex  shed  as  an  initial 
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condition  at  a  given  i)oint  on  the  particle  trajectorx’.  and  then  the  ua\f  p<uh(  t' 
can  be  superposed.  The  coin|)utatioiial  l>ox  should  he  long  enough  to  include  ihe 
effects  of  the  first  few  wave  packets  convected  downstream  while  the  next  lew  are 
incipient  near  the  particle. 

The  initial  conditions  due  to  vortex  shedding  are  periodic  in  time,  llowewi. 
we  cannot  simply  take  a  Laplace  transform  in  time  to  include  all  the  \•oHl((■^ 
because  the  particle  itself  is  moving  in  the  vertical  direction  '  the  dinxiion  of 
integration  of  the  stability  equation.  Hence  the  computations  will  be  carried  out 
in  the  time-marching  scheme  outlined  above. 

4.2.  Further  work 

Once  the  moving  particle  is  modeled  appropriatelv.  the  wave  packet  computations 
have  to  be  interfaced  with  the  trajectory  calculations.  The  end  result  sought  is 
an  efficient  algorithm  for  computing  the  disturbance  due  to  a  moving  particle 
of  any  size  and  specific  gravity.  Disturbances  due  to  two  or  more  particles  ma\ 
also  be  computed  and  the  various  wave  packets  may  be  sujrerposed  to  find  the 
resulting  disturbance.  If  necessary,  the  effect  of  nonspherical  partic  les  can  al.so  be 
studied.  Since  it  uses  a  linear  analysis,  this  computational  method  is  limited  in 
that  it  cannot  handle  nonlinear  interactions  between  the  particle  ami  the  bmimlary 
layer  flow.  To  study  such  interactions  of  strong  disturbances  it  may  be  necessarv 
to  use  a  direct  numerical  simulation  technique,  which  is  outside  the  scope  of 
this  project.  In  the  next  section  we  consider  the  stability  of  a  wake  in  a  two- 
dimensional  boundary  layer.  This  analysis,  and  its  extension  to  three-dimensivmal 
wakes,  will  provide  insight  into  the  selection  of  initial  conditions  and  the  inlluem  e 
of  the  particle's  wake  on  the  development  of  the  wave  packet  it  triggers. 

5.  Stability  of  a  Wake  in  a  Boundary  Layer 

Hall  (1967)  performed  an  experiment  to  determine  the  mechanism  of  transition 
due  to  bluff  bodies  suspended  in  a  boundarv’  layer.  Transition  aj;peared  to  depend 
on  the  stability  characteristics  of  the  wake  rather  than  on  the  stability  character¬ 
istics  of  the  boundary  layer.  For  spheres  at  different  heights  from  the  wall  the 
boundary  layer  responded  differently.  In  this  section  we  examine  how  the  location 
of  a  wake  in  a  boundary  layer  affects  tne  stability  of  the  entire  flow.  .‘\  spatial 
stability  analysis  is  used.  Three  parameters  that  characterize  the  wake  are  con¬ 
sidered:  the  wake  velocity  defect,  the  wake  width  and  the  wake  height  above  the 
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5.1.  Background 

The  earliest  stability  calculations  for  a  two-dimensional  wake  were  |)ei  lornied  i)\ 
Hollingdale  ( 1 940)  and  .McKoen  (1956).  They  used  analytical  methods  to  imd  tin 
eigenfunctions  for  neutral  disturbances  for  large  Reynolds  number  wakes.  Sato  ainl 
Kuriki  (1961)  reported  an  experimental  and  a  numerical  study  of  an  artiiicially- 
excited  transition  of  a  wake  behind  a  thin  flat  plate.  Based  on  their  ol:)servalions. 
they  classified  the  transition  2one  into  three  subregions: 

1.  a  linear  region,  where  the  two-dimensional,  sinusoidal  \  i-loc  iiy  IIik  t  ual  loii' 
amplify  exponentially  downstream,  and  different  freciuencies  aic  clainpcd  oi 
amplified  independent!}' 

2.  a  nonlinear  region,  where  the  grow’th  rates  deviate  from  Ireiiig  exponential 
due  to  nonlinear  effects,  and  different  spectral  components  ma\'  inieracl. 
and 

3.  a  three-dimensional  region,  where  the  disturbance  is  no  longi.n  periodic,  and 
spanwise  fluctuations  ultimately  cause  the  wake  to  become  turbulent. 

In  Sato  and  Kuriki's  numerical  approach,  a  parallel-flow  approximation  wiih  a 
Gaussian  mean  velocity  distribution  was  used  to  describe  the  wake.  The  Ha}  leigli 
equation  was  solved  for  temporally-amplified  disturbances,  to  obtain  the  eigenval¬ 
ues  Q  and  c.  To  compare  the  calculated  amplification  rates  oc,  with  experinu'n- 
tal  observations,  the  temporal  results  were  transformed  into  spatial  growth  rates 
acijcr  using  the  phase  velocity  of  the  waves.  Sato  and  Kuriki's  calculations  and 
experimental  data  agreed  well  in  the  linear  region.  This  confirms  that  a  linear 
stability  analysis  is  indeed  applicable  to  the  wake  when  the  disturbances  are  small. 

Mattingly  and  Criminale  (1972),  noted  that  a  spatial  stability  analvsis  gi\'es 
a  more  realistic  picture  of  flow  stability.  They  performed  a  stability  analvsis. 
in  which  both  frequency  and  wavenumber  are  complex,  to  study  the  near-wake 
region.  They  compared  their  calculations  with  experimental  observations  for  the 
wcke  of  a  thin  airfoil  in  a  water  tank.  They  studied  both  mode  I  (sinuous  mode) 
with  symmi,trica!  vertical  velocity  disturbances  and  mode  II  (varicose  mode)  with 
anti-symmetrical  ones.  They  found  that  mode  1  disturbances  are  predominant 
in  the  wake.  They  also  found  that  the  von  Karman  vortex  street  formed  Ischiml 
streamlined  bodies  has  its  origins  in  the  near-wake  region.  This  is  in  contrast  t(j 
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the  proposal  by  Sato  and  Kuriki  that  the  sortex  street  fonns  in  tne  non-nncni 
region  when  a  single  row  i  ;  spanwise  vortices  breaks  u])  into  two  rows  dWci)  iium 
the  centerline.  This  may  be  connected  to  the  absolute  instability  ol  wakes  in  tlic 
very  near- wake  region. 

The  terms  absolute  and  convective  instability  describe  the  impulse  resi)onsc 
of  a  flow.  If  an  impulsive  small  disturbance  grows  exponentially  at  the  source  lo¬ 
cation  in  time  the  flow  is  termed  absolutely  unstable.  For  a  convectively  unstabh' 
flow  the  transient  is  converted  away  from  the  source  and  leaves  the  flow  ullinialcl> 
undisturbed.  Huerre  and  Monkewitz  (1985)  demonstrated  the  distinction  between 
absolute  and  convective  instability  m  inviscid,  parallel  shear  flows  by  studying  the- 
branch-point  singularities  of  the  dispersion  relation  for  complex  fiet|uen(  and 
wavenumbers.  According  to  the  Briggs-Bers  criterion.  develo|>ed  in  the  >tud\  ol 
plasma  instabilities,  a  flow  is  convectively  unstable  when  the  branch  point  sin¬ 
gularities  are  in  the  lower  half  complex-frequency  plane.  Huerre  and  .Monkewitz 
(19SS)  also  numerically  verified  the  Chomaz.  Huerre  and  Redekopp  ( 1988)  iiukIcI 
for  the  sequence  of  transition  in  the  cylinder  wake.  The  critical  Reynolds  number 
for  transition  from  stability  to  convective  instability  (based  on  the  cylinder  diam¬ 
eter  and  freestream  velocity)  is  5,  and  the  first  local  absolute  instability  a))pears 
at  a  Reynolds  number  of  25.  As  the  Reynolds  number  is  further  increased  to  -IT  a 
global  temporal  instability  sets  in  and  leads  to  the  formation  of  the  Karman  vortex 
street  with  a  shedding  frequency  close  to  the  linear  global  response  frequency. 

In  the  present  study  we  consider  the  convective  stability  of  a  two-dimensional 
wake  embedded  in  a  boundary  layer.  The  strength  and  location  of  the  wake  are 
varied.  We  examine  whether  the  flow  stability  is  controlled  b}-  the  wake  or  ihe 
boundary  layer  as  a  function  of  changes  in  the  w'ake  properties. 

5.2.  Governing  Equations 

We  assume  the  basic  flow  to  be  tw'o-dimensional  and  a  function  of  y  only.  Tlie 
disturbance  velocities  u  and  v  and  kinematic  pressure  p  are  functions  of  r.  y  and  I . 
and  are  small  compared  to  the  basic  flow,  so  that  the  equations  of  motion  can  lx- 
linearized.  On  nondimensionalising  the  equations  (using  the  free.slrcam  wlociix 
(7(50  and  a  characteristic  length  such  as  the  wake  half-width), and  suiHracling  ilic 
mean  flow,  we  get  the  equations  governing  the  disturbances 

fir  +  Cy  =  0  1  18) 
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Following  the  conventional  normal  mode  approach  in  hydif-dynamif 
theory,  we  seek  complex  solutions  of  the  form 


where 


u{x,y.t)  = 

vix.y.i)  =  r(i/)c'<“"-"> 
pix.y.i)  = 

u  =  ^(f/  +  it')  V  =  ^(£’  +  v'}  p  =  -f  fr  j 


I  o  1  ) 
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and  the  superscript  *  denotes  complex  conjugate. 

Eliminating  the  pressure  term  and  using  the  continuity  tqualiun  leads  to  a 
single  differential  equation  for  v{y):  the  Orr-Sommerfeld  equation. 


-  q2)2  _  (C:  _  -  q2)  +  =  0 


a 


(■>'>} 


where  D  denotes  d/dy.  If  the  viscous  terms  are  neglected  we  obtain  the  Hayleigli 
equation 


(U  -  ~){Dh^  -  a^v]  ~  D^Uv  =  0 

Q 


m] 


with  boundary  conditions 


t;(0)  =  0,  Dv{y)-P  av(y)  -*  0,  as  y  — ♦  oc  (.■)7) 

w’here  a  =  +  ia,  is  the  wave  number  and  u;  is  the  frequencv  of  the  disturbance. 

For  a  two-dimensional  Gaussian  wake  in  a  flat-plate  boundary  layer,  the  mean 
velocity  distribution  is  the  superposition  of  the  two  profiles 

^(y)  =  -  u^.expi-a^iy  -  y^.)^)  (.jS) 

where  Ub{y]  is  the  Blasius  profile,  Uu.  is  the  wake  centerline  defect,  a^.  is  the  wake 
width  parameter  and  yxt,  is  the  wake  centerline  location. 
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5.3.  Method  of  Solution 

The  Rayleigh  equation  and  the  boundary  conditions  constitute  a  lioniogeneuu^ 
boundary  value  problem.  To  determine  a  solution  we  use  a  shooting  technicjue. 
The  mean  proiile  is  found  by  integrating  the  third-order  Blasius  e(|ualioii  as  well 
as  the  second-order  Rayleigh  equation  using  an  adaptive  step-size  Hunge  Kutta 
routine  [Press  et  al  (19S6)].  For  a  spatial  stability  calculation,  the  lrequc!ic\  _■  is 
real  and  the  wavenumber  a  is  complex. 

For  a  given  a,’  and  a  specified  set  of  wake  parameters,  the  system  of  li\  (‘  lit  si - 
order,  ordinary  differential  equations  is  integrated  from  the  wall  to  the  outer  e(lgc 
of  the  boundary  layer  for  two  guessed  values  of  a.  A  (complex!  secant  meiiiuil  |s 
used  to  obtain  subsequent  gues.ses  for  o  till  convergence  [Press  ei  al  (  lo 

obtain  a  reasonable  initial  guess  for  the  eigenvalue o.  we  use  t  he  results  of  Sat i;  and 
Kuriki  with  the  wake  outside  the  boundary  layer  (u^  =  0.7.  (i^  =  U.OOd.  i/„  =  (i.U) 
and  move  the  wake  down,  in  small  steps,  to  the  desired  wake  parameters.  1  ypn  al 
wake  parameters  used  in  this  studv  are  =  0.5,  a^,  =  3.24.1.  =  2.4  (Figure 

12). 

Once  eigenvalues  a  are  obtained  for  a  range  of  u.',  a  plot  of  the  spatial  growth 
rate  a;  versus  the  frequency  shows  the  most  unstable  frequency,  that  is,  the  Ire- 
quency  with  the  most  negative  o,.  For  this  mode,  eigenfunctions  e(.v)  «iid 
(from  continuity)  can  be  calculated.  The  energy  production  term  //|  can  bt*  com¬ 
puted  from  the  Reynolds  stress  term  —uv 

=  -uvVy  =  -^3?(ui-‘)Fv  (5‘)) 

This  term  is  a  measure  of  the  energy  transferred  from  the  mean  flow  lo  the  dis¬ 
turbances. 

5.4.  Results  and  Discussion 

To  study  the  effect  of  wake  parameters  on  the  stability  chaiactcristics  we  hist 
consider  Sato  and  Kuriki's  wake  added  to  the  boundary  layer  profile  far  Irom  the 
wall:  =  O.T.a^  =  0.7,  y^.  =  6.0.  .A.n  inviscid  analysis  is  performed  first. 

5.4.1.  Wake  Velocity  Defect 

For  a  fixed  a^.  and  as  the  wake  defect  is  decreased  from  0.7  to  0..5.  the 
rates  of  amplification  decrease  for  all  frequencies,  since,  in  the  limit  of  no  wake, 
the  Blasius  profile  is  stable  in  the  inviscid  analysis. 
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Figure  12:  Mean  velocity  profiles.  The  solid  line  shows  a  typical  wake  at  y^.  =  2.-1 
and  the  dashed  line  is  Sato  and  Kuriki’s  (1961)  wake  located  outside  the  boundai  x 
layer  at  y^,  =  6.0 
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Figure  13:  Variation  of  growth  rates  with  change  in  wake-widlli  in  h  l;ouiKlrtr.\ 
layer  with  =  0.5,  =  2.4 

5.4.2.  Wake  Width 

For  a  fixed  and  j/u,,  as  the  wake  width  parameter  is  increased  from  O.T 
to  3.8,  that  is,  as  the  wake  becomes  more  slender,  the  growth  rales  increase  for 
higher  frequencies  (Figure  13).  This  is  because  of  the  presence  of  sleeper  \elocily 
gradients  Uy  in  the  flow,  which  results  in  a  larger  positive  //j  in  the  viciniiv  of 
the  wake.  This  means  there  is  a  net  transfer  of  energy  from  the  mean  flow  to  the 
disturbance. 

5.4.3.  Wake  Height  in  the  Boundary  Layer 

For  a  fixed  u^,  and  Ou,,  as  the  wake  centerline  location  y^.  is  decreased  from  6.0  to 
2.4,  the  growth  rates  increase  for  lower  frequencies,  decrease  for  higher  freciiieiu  i<•^ 
and  so  the  most  unstable  frequency  decreases  (Figure  14).  .Now.  as  the  wake 
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Figure  14:  Variation  of  growth  rates  with  change  in  wake-height  in  a  Ijounclaiy 
layer  with  Uu,  =  0.5.  =  3.243 

approaches  the  wall,  the  “effective”  wake  defect  decreases  for  a  fixed  n,. .  since- 
the  local  mean  velocity  for  the  boundary  layer  is  lower  near  the  wall.  So  the 
inflectional  instability  decreases. 

5.4.4.  Viscous  Effects 

A  comparison  of  the  inviscid  Rayleigh  equation  solutions  with  the  solutions  to  thc 
Orr-Sommerfeld  equation,  we  find  that  the  addition  of  the  viscous  terms  docs  ikjI 
significantly  affect  the  stability  characteristics  of  the  wake.  This  means  that  ihc 
dominant  mechanism  of  instability  is  the  inviscid  mechanism  due  to  the  infieciional 
wake  profile.  The  details  of  the  boundary  layer  profile  near  the  wall  do  not  seem 
to  affect  the  stability  calculations.  Thus,  we  can  simply  put  the  wake  in  any  shear 
profile,  such  as  the  hyperbolic  tangent  profile,  having  the  same  slope  at  the  wake 
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Figure  15:  Variation  of  growth  rates  with  change  in  wake- velocity  deiecl  in  a 
hyperbolic  tangent  shear  layer  w’ith  a^.  =  3.243  and  /5  =  1.0  and  9.0 

centerline  as  the  boundary  layer  profile.  We  consider  the  Ijasic  flow  to  lie 

^^{y}  =  +  tanh{y))-  u^exp{-a^.(yf)  (GO; 

Here  the  parameter  for  the  location  of  the  wake  centerline  is  effectively  replaced 
by  the  parameter  0  for  the  slope  of  the  shear  profile  at  the  wake  centerline.  3'his 
has  the  added  advantage  of  giving,  in  the  two  limits  of  0  and  .  .he  profiles  stud¬ 
ied  by  Michalke  (1965)(0  =  1.0,  Uu,  =  0.0)  and  by  Sato  and  Kuriki  (1961  )(large  .1. 
and  Uu,  =  0.5).  Figure  15  shows  the  spatial  growth  rates  for  the  wake  in  a  shear 
layer  profile  with  a  wake  width  parameter  a^.  =  3.243  and  varying  and 
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5.5.  Future  Work 

The  problem  of  the  stabilit}'  of  a  wake  in  a  shear  flow  can  be  .solved  anai\  iical)\ 
using  a  perturbation  technique.  The  mean  flow  is  regarded  as  a  pei  lui  b.ii  um 
series  in  the  shear  parameter,  that  is.  the  shear  flow  is  treated  as  a  first -uidci 
correction  to  the  basic  wake  flow.  The  eigenvalue  (wavenumber)  is  also  e,\|;aiKletl 
in  a  perturbation  series.  When  these  expansions  are  substituted  into  the  iwc^- 
dimensional,  inviscid  equations  of  motion,  we  recover  the  Rayleigh  equation  loi 
the  unperturbed  wake  flow  as  the  zeroth-order  equation.  The  first-order  ecjuation 
is  a  non-homogeneous.  self-adjoint  equation  with  the  same  linear  ojrerator  a.s  lix' 
Rayleigh  equation.  This  can  be  solved,  using  the  soKability  condition,  to  obi.nn 
a  first-order  correction  to  the  eigenvalue  due  to  the  prescnc*'  of  slieai.  In  the 
two-dimensional  case,  the  results  may  be  compared  with  the  direct  nninerical 
calculations  reported  above.  This  technique  may  then  be  exteinled  to  the  tliier 
dimensional  problem  for  a  spherical  wake  in  a  two-dimensional  boundai  \  layei 
Barston(  1991 )  has  derived  criteria  for  the  instability  of  jrlane  paiallel  pio- 
files  with  multiple  inflection  points,  which  are  an  extensions  of  the  Ravleigh  and 
Fjertoft  criteria.  Hence  we  can  analytically  examine  the  sf'bilitx’  of  a  wake  in 
a  shear  flow  by  looking  at  the  relative  positions  of  its  two  inflection  points  in  a 
given  profile. 

5.6.  Conclusions 

The  stability  characteristics  of  a  wake  in  a  boundary  layer  have  been  detci  imned 
as  the  wake  approaches  a  wall.  These  calculations  help  provide  some  insight  into 
the  effect  of  particle-generated  wakes  in  laminar  boundary  layers.  The  results 
of  our  present  calculations  will  be  compared  with  ongoing  experiments  to  verily 
when  a  linear  or  weakly  nonlinear  stability  analysis  may  be  used  to  describe  such 
wake  flows. 
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